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GENERAL INTRODUCTION 

This document represents the Final Report on technical work 
accomplished under Contract NASV1-2953 for the period 1 duly 1976 to 
30 dune 1977. 

The organization of the Report is as follows. This General 
Introduction Section is followed by a brief section listing all of the 
oral presentations of work under the contract to scientific meetings, 
and of all technical papers either submitted to or accepted by technical 
journals. 

The technical work accomplished is then presented in five Appendices, 
each in the format of a technical paper, and therefore almost self- 
contained. Each Appendix contains its own summary and/or conclusions 
section. 

These Appendices are, in order, titled Oblique Whistler Instabilities,^ 
Current-Limited Electron Beam Injection, Three-Dimensional Ion Sound 
Turbulence, Theoretical Aspects of Sounder Antenna Operation, and Whistler 
Modes in Bow Shock Structures. - , 

The first three Appendices represent experimental work performed in 
TRW's LAMPS facility under the direction of Dr. Reiner L. Stenzel . The 
final two Appendices are theoretical work. The work on sounder antenna 
theory was done by Dr. R.W. Fredricks, while the work on whistler modes 
in bow shocks was performed by E.W. Greenstadt. All work reported here 
was performed under the supervision of Dr. R.W. Fredricks, the Principal 
Investigator. 
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ORAL PRESENTATIONS AND PUBLICATIONS 

The following oral presentations of work under this contract were 
made by the person whose name is underlined: 

1. R.W. Fredricks and R.L. Stenzel : “Observations of Return-Current- 

Driven Ion Acoustic Instabilities in a Beam-Plasma System," paper 
SM-34, Annual Fall Meeting, American Geophysical Union, San Francisco, 
December 1976. 

2. R.L. Stenzel and R.W. Fredricks: "Observation of Beam-Generated 

VLF Hiss in a Large Laboratory Plasma," paper SM-35, Annual Fall 
Meeting, American Geophysical Union, San Francisco, December 1976. 

3. R.W. Fredricks : "Laboratory Experiments of Interest to Space Plasma 

Physics," paper SM-49C, Annual Spring Meeting, American Geophysical 
Union, Washington, D.C., June 1977. 

The following technical papers have been accepted for publication 
in the journals listed: 

1. R.L. Stenzel : "Observation of Beam-Generated VLF Hiss in a Large 

Laboratory Plasma," accepted by Journal of Geophysical Research- 
Space Physics , 1977. 

2. R.L. Stenzel , "Experiments on Current-Driven Three-Dimensional Ion 
Sound Turbulence, Part I. Return-Current-Limited Electron Beam 
Injection," accepted by Physics of Fluids , 1977. 

3. R.L. Stenzel , "Experiments on Current- Driven Three-Dimensional Ion 


Sound Turbulence, Part II. Wave Dynamics," accepted by Physics of 
Fluids, 1977. 
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OBSERVATION OF BEAM GENERATED VLF HISS IN A LARGE LABORATORY PLASMA 

R.L. StsnzeV* 

TRW Systems, Redondo Beach, CA 90278 
ABSTRACT 

The predicted Instability of obliquely propagating whistler waves in 
a plasma penetrated by an electron beam Is verified in a laboratory ex- 
periment. The observations support the model for the generation of 
auroral hiss and compare favorably with ground and satellite obervations 
of VLF hiss. In contrast to the conventional smal 1 -diameter laboratory 
beam-plasma systems the present device is large compared to the characteristic 
whistler wavelength (\ - 2-4 cm, beam and plasma diameter 45 cm, length 
250 cm, resp.). Unstable whistlers can therefore, propagate and grow 
oblique to the beam over many wavelengths before encountering the 

plasma boundaries. When a cold (leV), energetic (40 eV) , low density 
8 “5 

(n^- 10 cm ■*) electron beam is injected along the magnetic field (B q =75 G) 
into an initially cold (kT e = 1/4 eV), dense (n e s 10^ cm background 
plasma the spontaneous emmission of broadband rf noise ranging from the 
lower hybrid frequency 01 ^ to the electron cyclotron frequency is observed. 
Performing narrowband two-dimensional cross-correlation measurements 
between two small antennas the rf noise is identified from its dispersion 
characteristics u(k) as whistler waves propagating near the oblique resonance 
angle, 9 - arccos (w/u c ) . The parallel phase velocity is close to the beam 
velocity (ta/kj j < v^) . Although the instability amplitude is far above the 
thermal noise level, it Is not large enough to give rise to significant 
nonlinear effects. The instability saturates when the higher frequency 
short wavelength electrostatic Instability near the electron plasma 
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frequency (u- »&£) forms a p-Jateau in the beam di stri but ion . Time and 
space resolved electron pitch-angle distribution function measurements 
are performed. The propagation of the whistlers out of the instability 
region and refraction into eleccromagneti c wh i stl ers is observed. Finally, 
in the unstable system, coherent whistler waves are launched from a small 
antenna, in the direction oppos i ng the beam the wave energy diverges along 
the resonance cone of angle 6 = arcs in Ww ) just as in a stable 

» C 

plasma. But j_n the direction of the beam, an unstable oblique whistler 
wave with a well-defined wave number is excited. The wave normal is 
at an angle 9 = + arccos (w/wj with respect to B Q , and the amplitude increases 
. aw . g y from the source at the complementary angle. Thus out of the broad k-vector 
spectrum excited by a point source the most unstable mode of the system 
predominates. These observations are relevant to the radiation characteristics 
of satellite antennas In regions of electron precipitation. 
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1. INTRODUCTION 

The phenomenon of auroral hiss has received considerable attention. 

The present model for the generation of VLF hiss is that of a collective 

interaction of an electron beam with a dense background plasma [Swift and 

Kan . 1975; Maqgs . 1976], An electrostatic instability arises when the beam 

velocity matches the parallel phase velocity of an oblique whistler wave 

[ Brinca , 1972] which can be very small when the wave propagates near the 

oblique resonance. Many observations from ground, rockets and satellites 

are consistent with the above model [He! 1 iwel 1 . 1965: Laaspere et al . . 1971; 

Gurnett and Frank . 197 2 ; James . 19733* The measurements usually consist 

of noise frequency spectra and electron pitch angle distribution functions. 

There are fundamental limitations of measuring wave normals or rays from 

a single observation point. These essential aspects of the instability are 

usually Inferred so that a conclusive proof of the mechanism cannot be given. 

it is suggested that a carefully designed laboratory experiment can supply 

this and other complementary information and thereby contribute to the 

understanding of the VLF hiss phenomenon. 

This paper d 'scribes a laboratory experiment on a beam-plasma instability 
under conditions for which the theory of VLF hiss is applicable. An 
energetic electron beam is injected into a cold, dense, almost collisionless 
magnetoplasma. Beam and plasma dimensions are so large that obi ioue 
whistier wave propagation and growth are readily observable. In this aspect 
the present experiment fundamentally differs from most previous laboratory 
experiments with small diameter beam-plasma systems which are restricted 
to one- dimensional wave propagation [ Malmberg and Wharton. 19^9; Bol 1 inger 
et al . . 1974; Gentle and Lohr . 19731* The present situation models in many 
aspects the conditions inside an auroral arc. It does not apply to the 
wave propagation and refraction effects outside the arc since beam and 
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plasma diameters are identical. Whistler wave propagation in the absence 
of the beam have been reported earlier [ Stenzel . 1976 a, b, c] . 

The observations show that the beam generates broadband rf noise 
(ui^<isi<w c ) . By analyzing the frequency and wavenumber spectrum, i.e. the 
dispersion oj(k), the noise is identified as whistler waves propagating 
near the oblique resonance cone with parallel phase veloci ty si ightly less than 
beam velocity. This is a direct confirmation of the essential prediction 
from the instability theory. 

The dynamics of the beam and background plasma are investigated. It 

is concluded that the simultaneously occuring short wavelength electrostatic 

instability near the electron plasma frequency (oj >« )is responsible for 

p c 

the scattering of the beam and the saturation of the oblique whistler ins- 
tab 1 1 i ty . 

Finally, wave propagation experiments are performed by launching 

coherent “test 11 waves from a small antenna in the unstable region. In the 

direction opposing the beam the point source excites at a given frequency 

a broad wave number spectrum, k (0 ) , which interferes to give the familiar 

resonance cone pattern [ Fisher and Gould . 1971, Kuehl , 1962]. But in the 

direction wi th the beam the most unstable mode with kj j > w/v^ and 

kj/kj | = tan6 = (a) 2 /u 2 -1)^^ dominates, i.e. a single oblique whistler mode 

is excited whose wave normal angle is close to that for oblique resonance 

and whose amplitude grows in the ray direction at nearly the complementary 
angle. 

The connection and relevance of the laboratory results on observations 
in space are discussed in the last section. 


.... -| - r 
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!!. EXPERIMENTAL SET-UP AND MEASUREMENT TECHNIQUES 




The experiment is performed in a laboratory plasma source designed 

to produce a large diameter uniform collisionless magnetoplasma. The 

device which has previously been described in more detail [Stenzel , 

1976 a, c] consists briefly of a dc discharge between a solid oxide-coated 

cathode of approxim. 1/2 m diam. and adjacent grid anode [Fig. 1(a)]. The 

plasma is confined by a uniform axial magnetic field (B = 75 G) and 

terminated after 2.5 m by a floating end wall. As shown in Fig. 1(b), 

the discharge is pulsed repetitively. The first pulse produces an after- 

9 11 ”3 

glow background plasma which is cold (kT = kT, = 0.25 eV), dense (Ur <n <10 cm , 

ti | C 

-k 

1 < to /tu < 10), weakly collisional (electron-neutral collisions v /to- 10 , 

pc' en c 


cou 


_ *2 

lomb collisions v ./ui = 10 J ) and free of low or high frequency ins- 

c 


tabilities. At a desired afterglow time a second voltage pulse is applied 
to the cathode whereby a beam of fast electrons is injected along B q into 
the afterglow plasma. Beam and background plasma parameters are determined 
from plane Langmuir probes which can be positioned radially and axially. 
Independently, the density is determined from dispersion measurements of 
whistler waves, Langmuir waves and Bernstein waves, and from resonance 
cone measurements. High frequency signals are detected and/or excited 
with two identical short wire rf probes (5 mm long exposed center wire 
of 50 fi semirigid coax cables) one of which is movable along and across 
B . Due to the pulsed nature of the experiment ail data are sampled and 
averaged over many highly reproducible events. By varying the sampling 
time and probe position time and space resolved data are obtained so that 
one can follow the evolution of the beam plasma Interaction from turn-on 
to steady state. 




T 


I 
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111, PLASMA AND BEAM PROPERTIES 

First, the background plasma properties are presented. Fig. 2 shows 
the density decay of the afterglow plasma. The density is measured by 
different independent techniques which, In general, show good agreement. 

The plasma decay is due toambipolar diffusion where the dominant loss is 
due to surface recombination on the cathode, anode and axial end wall. 
Radially, the plasma is separated from the chamber wall by a 20 cm wide 
gap in which the density is negligibly small. 

The electron temperature during the active discharge is approximately 
kT g = 2eV. in the afterglow, at 2 10 msec, the temperature has decreased 
kT = kT. - 0.25 eV. The electron temperature is determined from 

£ I 

Langmuir probe measurements and, independently, from the dispersion of ion 
acoustic waves. Landau damping of ion sound waves yields the estimate of 
the ion temperature. Other important parameters of the background plasma" 
are 1 isted in Table 1 . 

The electron beam Is injected into the background plasma whose density 
is determined by the afterglow time t . Beam and background plasma strongly 
interact and change their initial properties to new steady-state values. 
Figure 3 shows the electron distribution function sampled at different 
times t g after the step-like turn on of the beam voltage at an after- 
glow time t = 15 msec. The distribution function is obtained by dif- 

3 

ferentiating the current-voltage characteristics [f e (V) “ dl/dV] of a 
one-sided plane Langmuir probe located at a distance z = 7 cm from the 
anode with its surface normal pointing toward the beam source. One can 
see that initially (0 < < 30 ysec) a cold electron beam is injected 

into a cold background plasma (kT^ < kT^ « eV^) . This system is- highly 
unstable to electrostatic instabilities. With increasing time the beam 
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broadens (t g - 40 psec), forms a plateau (t g = 50 psec), and finally 
(t g s 60 psec) becomes the tail of a distribution with continuously 
decreasing slope- 3f /3v < 0. On the same time scale, the background 

electrons are strongly heated such that both distributions finally merge. 

One also observes that beam electrons are accelerated to energies larger 
than their initial injection energy. This fact as well as the short 
penetration length of the beam compared to the electron mean free path 
indicate that the beam is scattered by electrostatic instabilities rather 
than by collisional processes. 

On a longer time scale, however, the beam electrons undergo collisions 

with the low pressure background gas. When the beam energy exceeds the 

ionization potential (E-. 11 14 eV in krypton), the plasma density increases. 

Ionization phenomena become noticeable for t > 50psec, When a 40 eV beam 

s 

pulse Is applied for approx. 500 psec, which are the typical conditions of 

1 1 -3 

the discharge pulse, the density rises to the value at t =0 (n =2 x 10 cm / . 
For the purpose of the present experiment, however, the initial period 
(0 < t < 50 psec) of the beam injection is of particular interest since 
unstable distribution functions are clearly observable and collective 
processes dominate. This period Is long compared to the time scale of 
the rf i ns tab i 1 ities. 

The spatial properties of beam and plasma distributions at a given 
time t = 10 psec are shown in Fig. 4. With increasing distance z from 
the injection point the beam broadens and the background electron tem- 
perature and density build up. Unstable beam distributions only persist 
up to distances z £ 20 cm from the cathode. The penetration depth depends 
on beam speed and density but not on neutral pressure. 
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Typical radial and axial density profiles in the vicinity of the cathode 
where the oblique whistler instability is investigated are shown in Fig. 5- 
The density variations appear small enough to consider the medium as 
uniform. The magnetic field is spatially uniform to Within +_Q.$% and 
temporally constant to within 0,1%, 

The probe measurements presented so far give only information about the 

parallel electron distribution function. By rotating the one-sided plane 

Langmuir probe as sketched in Fig. 8, the pitch angle distribution can 

be estimated. Figure 6(a) shows the distribution function at different 

angles <f> between the surface normal and the magnetic field. The measurements 

are taken near the cathode (z = 7 cm). At ej> = 0 the probe faces the beam 

and a large beam current is collected when jV^lc^. With increasing angle 

<f> the peak of the beam distribution shifts to smaller retarding potentials 

! V I and the beam current decreases until no beam is detected when the 
1 P l 

probe surface is parallel to B q . This behavior identifies the beam as 
being cold and parallel to B q . The beam temperature kT^ can be estimated 
from the energy spread AV^ near onset, kT^ - {l/4)eV^ (AV^/V^) s 0.5 eV. 

When the probe is inclined at an angle <J>, the normal beam velocity component 
is v b costfi and beam particles are detected at probe voltages | | < cos <j>. 

The collected current decreases as cos [see Fig. 6(b)] because both the 
normal beam flux and the probe area projected along B q decrease as cosfjt. 

The background plasma, how- ver, exhibits a relatively isotropic distribution. 

Figure 7 shows a similar measurement as Fig. 6, but it is taken in the 
center of the plasma device (z - 125 cm) where the beam has been modified by 
electrostatic instabilities. The probe current increases exponentially in- 
dicating that the beam has been thermal ized ( kT^ - 8 eV) . The angular variation 

shows a continous decrease in energetic particle current with increasing angle <!> 
but little change in the energy spread. 
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F or <p = 90 ° the flux at a given voltage is more than an order of magnitude 
smaller than for <j> = 0, and for <j> = 180°, particles at the incident beam 
velocity cannot be detected within experimental accuracy. These observations 
Indicate that the beam is not scattered isotropi cal ly but mai nly . broadened 
along B q in the direction of injection. Some angular spread occurs since a 
small energetic perpendicular flux is observed. The distribution function 
is stable against high frequency electrostatic instabilities and exhibits 
no perpendicular temperature anisotropy. 

While Fig. 7 shows only the energetic particle distribution. Fig. 8 
Includes the behavior of the background plasma. The bulk electron temperature 
is approximately kT^ = 3 eV. However, due to the thermalized beam, the parallel 
distribution function is asymmetric. Since the electron saturation current 
is identical in both directions, the plasma is current free which is consistent 
with the fact that the discharge current flows only in the 1 cm wide gap 
between cathode and anode. The injected beam current passing through the mesh 
anode ( 70 % transparent) is exactly balanced by a return current to the anode. 
This current consists of background electrons drifting toward the anode. The 
drift velocity is smaller than the beam velocity by the ratio of beam to 
background electron density (n, /n =10 ). This small drift (m Vj /2 < lev) 

which is not even noticeable in Fig. 8 has no direct effect on the high 
frequency instabilities. However, since v d » c g = (kT^/mj) ' ion acoustic 
waves are driven unstable [ Mi khaijovskl , 197 ^] and the resultant strong 
turbulence [ Kadomtsev, 1965] scatters high frequency waves. The observations 
on ion sound turbulence driven by the return current are beyond the scope of 
the present paper and will be discussed elsewhere [Stenzel , 1977]* 
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IV. INSTABILITY ANALYSIS 

The onset of an Instability upon beam injection becomes apparent 

through the large rf noise level detected with electric or magnetic antennas 

inside the plasma. Fig. 9(a) shows the coincidence between rf noise and 

electron beam. The top trace is the rf potential fluctuation which, for 

purpose of display, has beeen heterodyned with a local oscillator at 

f = 100 MHz (ui /(u = 0 .A 7 ). The bottom trace shows the beam current pulse 
c 

applied at t = 10 msec in the afterglow of the main discharge. After an 
2 

initial trans i ent, the current is constant during the first 50 ysec and 
then rises as the density builds up. As outlined earlier, the beam current 
is balancec by the return current to the anode which is limited by the 
electron saturation current of the background plasma. The comparison bet- 
ween the two traces clearly shows that the rf is generated only in the 
presence of the electron beam. This coincidence holds for a wide range of 
frequencies extending from the regime of ion acoustic waves to that of 
electron plasma waves. 

The high frequency signals are applied to a spectrum analyzer with fast 
output response (300 kHz bandwidth). Figure 9(b) shows the frequency 
spectrum in the whistler wave regime sampled at t = 100 psec after turn-on 
of the beam pulse. The top trace is the rf potential from a short-wire 
antenna near the region of beam injection, the bottom trace is the magnetic 
loop signal detected in the middle of the device at z = 125 cm from the 
cathode. The broadband character of the instability Is apparent. The rf 
magnetic field spectrum is cut off beyond - O .85 which appears to be the 

result of refraction and damping of whistlers propagating from the beam 
injection region to the center of the device where the distribution is stable. 
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Although, in general, the shape of the whistler wave spectrum depends 

on beam and plasma parameters there is no particular structure noticeable 

-2 -2 -I -1 

near the characteristic frequencies, the lower hybrid (w^ - to . + w ce “ C |J 

= 1 MHz) and the electron cyclotron frequency (tu /2 tt = 210 MHz). 

! h ce 

Figure 10 shows a comparison of frequency spectra at two widely different 

densities [tD 2 /^ 2 = 2.8 in (a), to 2 /u 2 = 36 in (b) ] . In the low density - 
pc P c 

case, the frequency scan is extended beyond the whistler wave regime to 
show the presence of strong high frequency electrostatic instabilities near 
the electron plasma frequency. 

tn order to Identify the rf noise signals as waves, the spatial 

properties have been investigated with correlation techniques. Two 

identical movable antennas are used whose noise signals are passed through 

narrowband amplifiers, multiplied, time averaged on the high frequency time 

scale and then sampled and averaged on the longer time scale of the beam 

pulse. A typical cross correlation signal is plotted vs. axial position 

in Fig. 11(a). One can see that the noise correlates for a few wavelengths 

along the beam direction. By measuring the wavelengths for different 

frequency components, the dispersion u(kjj) is found. Figure 11(b) shows 

that the parallel phase velocity is slightly below but close to the beam 

velocity. This is characteristic for a Cherenkov type beam-plasma instability. 

The direction of wave propagation is determined from the measurement 

shown in Fig. 12. A delay line is inserted in one channel of the correlator 

[Fig, 12(a)] and the spatial correlation function is plotted at different 

delay times t [Fig, 12(b)], From the time-position diagram of stationary 

phase points of the correlation signal [Fig, 12(c)] one concludes that the 

waves propagate in the direction of the electron beam at a slightly smaller 
phase velocity than the beam velocity. 
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So far, only the component kj | of the wave vector has been determined. 
The perpendicular wave number is found by performing two-dimensional 
cross correlation measurements. Figure 13 (a) shows a map of the location 
of minima and maxima of the correlation signal in the r-z plane. One 
finds that the phase fronts are inclined with respect to the magnetic 
field. The wavenormal or k makes an angle 0 with respect to which is 
close to but below the angle for oblique resonance, 0 = arccos (co/w^) . 
According to cold plasma theory [ Bei 1 iwel 1 , 1965] the dispersion for 
oblique whistlers is given by 


.2 2 

k c 


to 


1 + 


0 ) 


to (to COS 0-to) 
c 


With the experimental parameters to, to , to given in Fig. 13 and the 

P c 

measured oblique wavelength A - 2.2 cm the calculated angle 0 is found 
to be 6 - 38.9° which is in good agreement with the measured value 0=38.7°. 

The angle for oblique resonance is 0 = ItO.S 0 . 

By measuring various phase front maps at different frequencies one 
finds that the propagation angle changes with frequency as shown in Fig. 13(b). 
Over a wide range of frequencies the rf noise Is thus identified as oblique 
whistler waves propagating near the oblique phase velocity resonance cone. 

Two further observations are worth pointing out. First, the sign of 6 is 
found to vary arbitrarily. At some frequencies the waves propagate from 
the interior of the plasma outward, at other frequencies they propagate 
inward. Thus, there is no constraint on the direction of wave propagation 
which is characteristic for an unbounded medium. Second, the amplitude 
pattern of the cross correlation signal shows a systematic frequency depen- 
dence. At lower frequencies (to/<u < 0.5) the correlation length along 

C 

B is larger than across B , at higher frequencies (w/w -> 1) the opposite is 
0 o c 

observed. This behavior is consistent with the fact that group and phase 
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velocity vectors are essentially perpendicular to one another near the 
resonance cone. The correlation amplitude is large within the group velocity 
cone which narrows toward lower frequencies. 

A measurement of the absolute value of the instability amplitude 
involves the knowledge of the coupling properties between antenna and 
wave. These are difficult to determine theoretically but can be estimated 
experimentally. For example, as described in the next section, a resonance 
cone is excited by applying a known small amplitude rf signal to one of 
the two identical short wire antennas. With the second movable antenna 
the absolute received signal along the resonance cone is measured. The 
coupling loss between the two antennas, extrapolated to vanishing probe 
separation, can be considered as twice the wave-antenna coupling factor 
for a single antenna (typical ly-22dB) . Similar techniques have previously 
been used for absolute measurements of electron plasma waves [ Kalmberg and 
Wharton , 1 967] • With this calibration procedure the absolute rf potential 
fluctuations are determined to be U r ^| “ 8 mV or e| <j> r ,p|/kT e s 5 x 10 
At these amplitudes nonlinear wave effects are not considered important. 

A second independent measurement of the wave amplitude has been 
performed in the middle of the plasma device. The oblique whistlers from 
the instability region near the cathode refract into nearly transverse 
whistlers which propagate along the entire plasma device and can be observed 
with electric and magnetic dipole antennas. The absolute wave magnetic 
field is measured with a shielded loop antenna whose response has been 
calibrated [ Stenzet , 1976 c] , Typical magnetic field strengths 10 ^ 

Gauss or |B r ^|/B Q - lo""*’ are observed which again shows that the wave remains 
in the linear regime. 
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V. TEST WAVE PROPAGATION 

The propagation of small amplitude phase coherent waves from antennas 
does not only provide a conclusive test whether the rf noise is due to in- 
coherent Cherenkov radiation or due to a beam-plasma instability but also 
provides further information on group and phase velocity directions, growth 
rates and saturation levels. 

When an rf siqnal (us... < us < us < us ) is applied to a small antenna in 
3 In c p 

a Maxwellian plasma, the radiation pattern assumes the shape of two symmetric 

cones. The origion is at the point source, the axis is long B q and the half 

angle is the group velocity resonance cone angle 6 - arcsin (w/us ) . This 

c c 

predicted pattern [Kuehl , 1 962.] has been observed both in the near zone 
[Fisher and Could , 1971] and in the far-zone [ Boswell and Gonfalone , 1975; 
Stenzel , 1976], The field distribution results from the interference of 
the wide range of oblique eigenmodes excited by an oscillating point 

source. 

In the beam-plasma system the radiation pattern shown in Figure 14 is 
found to be asymmetric along B q , In the direction against the beam a 
resonance cone is excited just as In the absence of the beam. Figure 
14(a) shows a single localized field maximum when the receiver probe is 
scanned axially at a fixed radial distance from the cone axis. A two- 
dimensional scan in the r-z plane [Fig. 14 (b)] shows that the field maxima 
diverge along two straight lines which, due to azimuthal symmetry represent 
a cut through a conical radiation pattern. The observed cone angle 8 c - 45° 
agrees with the expected value 0 = arcsin (to/to ) = 4y.2°. Since the wave 

energy flows along the cone through an increasing cross section, the wave 
amplitude rapidly decreases with distance from the exciter. Contour maps 
of equal relative rf potential | ^ (r, z) 1 are shown in Fig. 14(e). 


I 
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in the direction with the beam the interferometer trace, Fig. 14(b), 

reveals a well defined wave. Its amplitude grows with increasing distance 

8 

from the exciter, its parallel phase velocity (fXjj ” 3.7 x 10 cm/sec) 
is close the beam velocity. The phase fronts shown in Fig, 14(d) are 
oblique with respect to B"*" and the wave normal makes an angle 0=90-6 - + 45° 
with respect to B q . Thus, just as Is the spontaneous noise emission, the 
test wave Is an oblique whistler wave propagating near the resonance cone. 
The amplitude distribution shown in Fig. 14(f) indicates that the direction 
of maximum growth, i.e. the group velocity direction Is also oblique to 
and is nearly perpendicular to the phase velocity direction. These 
observations directly confirm the predicted oblique whistler instability 
in a beam-plasma system. 

The test waves are propagated tn the presence of the background noise. 
Since phase coherent detection is employed and test wave amplitudes above 
the noise level are chosen, the excited waves are readily distinguished from 
the background. The fact that the test wave grows indicates that the in- 
stability is not yet saturated which is consistent with the relatively small 
noise amplitudes measured previously. Although the Instability shows a 
broad spectrum of incoherent waves, it cannot be considered as a turbulent 
system in the sense that strong wave-wave or wave-particle interactions 
would take place. For the test wave. Fig. 14£f] gives an oblique e-fo]ding 
growth length k.^ s 3*5 cm and Fig. 14(d) gives the oblique wavelength 
X= Jbir/k^l.S cm. Because of the conical field distribution the measured 
normalized growth rate, k./k r = 8 x 10 , is smaller than that for a plane 

k 

wave so that a comparison with theory cannot be easily given. Nevertheless, 
one can estimate that a whistler at 6=0 =45° grows by at least 8 e-folding 

v 
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lengths ('v>70dB) after propagating obliquely through half the column of 20cm 
radius. Further convective growth is limited by the fact that the higher 
frequency instabilities (to— 01 ^} nonljnearly modify the beam. 

A continous transition from a resonance cone pattern to a growing wave 
due to a variation in beam voltage is shown in Fig. 15. The probe arrangement 
is the same as in Fig. 14(b). At low beam voltages there is a mixed pattern 
between cones and waves. The wavelengths can be as short as the half-width 
of the unperturbed cone; the refractive index for these oblique whistlers is 
n 2 150 whereas n - 4 for a parallel propagating whistler. Warm plasma 
modes [Fisher and Gould , 1971; Kuehl , 1973] are not observed. As the beam 
voltage is increased, the cone structure vanishes, the parallel wavelength 
and the growth rate increase. The latter may be due to an increase in the 
ratio of beam to background density with beam voltage. 

The importance of the simultaneously occurring electrostatic Instabilities 
near the plasma frequency has been stressed and their presence indicated in 
Figure 30(a). Figure 16 shows the behavior of a high frequency test wave 
near the most unstable upper frequency component of the noise spectrum. 

Figure 16(a) indicates that for a constant applied signal to the antenna 
large amplitude waves are only excited in the presence and direction of the 
beam. Although the waves are excited from a point source, the wave fronts 
shown in Fig. 16(b) are essentially plane and propagation is parallel to the 
beam. Due to the long interaction length with the beam, these waves exchange 
energy with the beam most efficiently, which leads to beam broadening and 
quenching of the instability. The temporal growth rate for the hydro- 
dynamic beam plasma instability is given by [ Mikhai lovski i , 1974] 
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Spatial and temporal growth rates are related by k. = y/v where the group 

^ 9 

-2 

velocity Is approximately given by the beam velocity. For n^/n^^ = 10 
one finds k.^ =1 cm which is much shorter than the whistler wave growth 
length component al'ong B^. 

The hydrodynamic instability corresponds to the case of a cold beam 
with a, /v, < (n./n ) where a, , v, are thermal and directed beam velo- 
cities, respectively. For a typical beam density n^/n^ ” ^ the beam 
energy spread is restricted to AV^/V^ = ^a^/v^ < 0.43. Figures 3 and 4 
show that this condition is satisfied in time up to ^ 40 psec after beam 
injection and in space up to ^ 15 cm from the beam source. Due to the 
injection of a cold beam the largest amplitude mode near f - 350 MHz in 
the high frequency spectrum of Fig. 10 (a) can be associated with the 
plasma frequency rather than the upper hybrid mode. The latter would 
imply highly oblique propagation which, according to Fig. 16(a) is not 
observed. The structures In the frequency range f c < f < f in Fig. 10(a) 
are not associated with characteristic frequencies of the plasma and are 
simply unstable modes on the beam with real frequency given by w - kv^. 

VI. DISCUSSION 

The laboratory observations show many similarities with phenomena 
observed on natural VLF hiss In space. The baste agreement is that VLF 
hiss is associated with the presence of an electron beam passing through 
a background plasma. A wide range of beam velocities gives rise to the 
instability. A broad band of frequencies is excited which shows little 
structure when the observer is inside the beam plasma system. The frequency 
spectra observed by James [1973] on lsis-1, which passed through regions 


of electron precipitation, are very similar to the spectra measured In the 
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laboratory. The various characteristic shapes of the frequency-time 
plots such as saucers and inverted v's shown by Gurnett and Frank [1972] 
are due to propagation effects between the instability region and the 
moving observer outside the arc. Likewise, observations from ground are 
strongly affected by wave reflection and refraction processes. The labora- 
tory experiment cannot model such propagation effects involving magnetic 
field and density gradients and long propagation distances. It does, 
however, provide much more inside into the generation mechanism of VLF hiss 
than spacecraft or ground observations. The theoretically predicted 
coherent Cherenkov-type instability of oblique whistlers is clearly iden- 
tified through direct measurements of wavenumber spectra, group and phase 
velocities, growth rates, instability amplitudes and electron distributions. 

The observations also show that the Cherenkov type Instability is dominant 
over other beam-plasma instabilities, for example the cyclotron-type 
instability associated with the anomalous Doppler effect, m=kj | v^-w c , 

[Mikhai lovski , 197^3 . No evidence is found for this instability, pre- 
sumably due to its smaller growth rate. Various authors [ Hashimoto 
and Matsumoto, 1976; Matsumoto and Yasuda, 1976; Cuperman et ai . , 1976] 
consider only wave propagation parallel to b£ which is too restrictive 
for an unbounded medium. 

The present experiment also provides insight into the evolution of 
the beam-plasma system. The electrostatic Instabilities of a cold beam 
plasma system are so strong that the beam rapidly spreads in velocity 
space. Host observations of the electron pitch-angle distribution In 
auroral arcs [ Reasoner and Chappell , 1973; Frank and flckerson , 1971] show 
beams with significant energy spread. Although the mechanism for beam 
generation Is not fully understood, it is possible that the energy spread 
is the result of electrostatic Instabilities as observed Tn the laboratory. 
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The finite spatial extent of the beam limits the convective amplitude 
growth of the oblique whistlers both in the laboratory and in space. 

However, wave propagation near the resonance cone is slow and the wave 
amplitude can grow well ab c'/e the level of incoherent Cherenkov radiation. 
For very low energy electron beams when the wavelength is within the 
resonance cone width, the instability in the laboratory vanishes. In space, 
Gurnett and Frank [1972] also observe a threshold electron energy level 
below which VLF hiss is not observed. 

The test wave propagation experiments which give useful information on 
the instability in the laboratory could also be performed in space. James 
[1975] has observed propagation near the resonance cone in i ntersateil i te 
radio transmissions between two of the Alouette~l sis spacecraft* If such 
"mother-daughter" type wave propagation experiments are performed in 
regions of electron precipitation, the transmission would be non- 
reciprocal exhibiting one-way amplification. 

VII* SUMMARY AND CONCLUSION 

The present paper describes an experimental investigation of a beam- 
plasma instability which is assumed to be the mechanism for generating 

auroral VLF hiss. The experiment is performed in a laboratory plasma 

of characteristic dimensions large compared to the growth length of the 

instability. The plasma parameters are appropriate for modeling the 

electrostatic instability of an oblique whistler wave driven resonantly 

by a parallel electron beam. The observations of frequency and wavenumber 

spectra, growth rate and instability level are consistent with theory and 

satellite observations in regions of electron precipitation. Due to the 

variety of diagnostic tools, a scaled laboratory experiment provides 

information not available from observations in space, and may thereby 

contribute to the understanding of the physics of VLF hiss. 
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FIGURE CAPTIONS 

Fig. I (a) Schematic diagram of plasma device. 

(b) Typical time sequence during one period {25 msec) of the 
pulsed beam injection experiment. 

Fig. 2 Density decay of the afterglow plasma measured by independent 
techniques. B q = 75 Gauss. 

Fig. 3 Electron distribution function I1 B q sampled at different times t after 
the turn-on of the beam pulse. The figure demonstrates the bro- 
adening of the beam and heating of the background plasma in 
time at a fixed distance z from the anode, t = 15 msec. 

fig. 4 Electron distribution function ll^ D sampled at different distances z 
from the cathode for a fixed Use t g after beam injection. The 
figure shows that beam and background heating occur over distances 
short compared to the electron mean free path (l > 90 cm), t =24 msec. 

£ 3 

Fig. 5 {a) Radial density profile of beam and background plasma-, 

n, = 5 x 10 7 cm ^ , n - 2 x 10^ cm z - 7 cm. t « 24 msec, 
b 3 o ’ 3 a 3 

B o “ 75 G ‘ 

(b) Axial density profile with and without beam, r - 0, t =24 msec. 
The figure shows that the beam-plasma system is reasonably 
uniform in the beam injection region where the instability is 
investigated. 

Fig. 6 Differentiated current-voltage characteristics of a plane one- 

sided Langmuir probe at different angles between surface 
normal and magnetic field (see sketch in Fig. 8). Beam 
retardation voltage and beam current tsee part {b)3 vanish 
as $ * 90° whereas the background electron current is rather 
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insensitive to probe orientation. The electron pitch angle 

distribution in the beam injection region (z = 7 cm) is that 

of an energetic, cold electron beam passing along through 

an isotropic background plasma, t = 24 msec, B q = 75 0* 

Fig. 7 Plane, one-sided Langmuir probe characteristics in the electron 

retardation region, log 1 vs. V , for different probe angles <j> 

P P 

at a large distance from the cathode (z = 125 cm). The beam 
has been mainly broadened in the direction of beam injection 
(kT^ = 8 eV) although some side scattering is observable. The 
beam shows no temperature anisotropy. 

Fig. 8 Langmuir probe traces in the middle of the plasma device (z = 125 cm) 

for the probe facing the beam (if» = 0) and opposing it (<j) = 180°), 

The plasma is free of a net current implying that the beam current 

is balanced by a return current consisting of slowly drifting 

background electrons (c g « Vj< v g ) . 

Fig. 9 (a) Rf noise at w/w =0.47 (top trace) and beam current pulse 

(bottom trace) vs. time, showing the cause-effect relationship 

between beam and rf emission, t =15 msec. 

a 

(b) Frequency spectrum of the spontaneous emission sampled at 

t = 100 ysec after beam injection. <j> f is the electrostatic 
s rT 

signal detected In the instability region (e | <f> ^ | /kT e = 5x10 5 ) 

and is the wave magnetic field of whistlers detected 100 cm 

away from the excitation region. 

I B /8 = 1Q“ 5 , S = 75 G, w/fc.- - 5. 

1 -rt 1 o o pc 

Fig. 10 Comparison of the instability frequency spectra at low densities 

[cu 2 /m 2 - 2.8 in (a)] and high densities [w 2 /u 2 = 36 In (b)], 
p p 3 p c 

The high frequency instabilities (m <«< m ) are also shown in (a). 

c p 

The observed spectra have the frequency response of the antenna folded in. 



29697-6002-RU-00 
Page A-26 

Fig. 11 (a) Cross correlation signal between two narrowband noise 

signals obtained from a fixed reference antenna and an 

identical antenna movable along . The parallel wavelength 

(A j | = 1.75cm) is much shorter than that of a whistler propagating 

along B (A = ?.9cm). di/to =0.9, u /to ~ 6, z ~ 10cm. 
o w c pc 

(b) Dispersion to vs. k j j of the spontaneous noise emission. 

The parallel phase velocity to/kj | is close to, but below the 
beam velocity. 

Fig. 12 (a) Schematic arrangement to measure the direction of wave 

propagation of the noise emission. 

(b) Correlation signals vs. position along the beam at different 
delay times t. 

(c) Time — position diagram for stationary phase points (circles = 
maxima, squares = rr.’nima) of the correlation signal. The figure 
shows that the parallel phase velocity points in the direction of 

i /n 

the beam with magnitude |v j < (2eV./m )* • 

p D € 

Fig. 13 (a) Phase fronts of the spontaneously excited waves in the 

r-z plane obtained from two-dimensional narrowband cross correlation 

measurements. The wave normal makes an angle 0 with respect to 

which is close to the oblique phase velocity resonance cone 

angle, 0 = arc cos (tu/oi ). 

c 

(b) Propagation angle 0 at different frequency components 

tu/ti) . The measurement shows that the broadband noise consists 
c 

of whistlers near the oblique resonance. ;= 50 V, t g = 30psec, 
z - 8cm. 
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14 Properties of small amplitude test waves excited from a point 

source In a beam-plasma system. An axial interferometer trace 

(a,b) shows a highly localized field for propagation against 

the beam but growing waves for propagation with the beam. A 

two-dimensional phase front map shows a resonance cone pattern 

against the beam (c) and oblique waves with the beam (d) . A 

contour map of constant rf amplitudes (arbitrary units) shows 

amplitude decay along the diverging cone opposing the beam (e) 

but oblique amplitude growth with the beam (f ) . Note that for 

the unstable mode u>/ki i S v. and v I v . Parameters: w /m =3, 

| | b g J- p pc’ 

t = 10 psec, t = 20 msec, V, - 38V. 

S 3 D 

15 Axial test wave interferometer traces at a fixed radial 

distance Ar from the exciter antenna at different beam 

energies eV^. A continuous transition from a resonance 

cone pattern (V^- 0) to a growing single whistler mode is 

observed, t - 2k msec, t = 30ysec. 
d 5 

16 Test wave propagation near the electron plasma frequency. 

(a) Interferometer signal ( 2 :) cos(kjjz) vs. z, showing 
that strong waves are only excited In the presence and 
direction of the beam, (b) Two-dimensional phase front 
map (maxima of interferometer signals). Although the 
exciter is a point source plane waves along the beam are 
excited which are the most unstable modes of the beam-plasma 
system leading to beam and background plasma heating. 
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Figure 11 
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Figure 12 
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Figure 14 
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Figure 15' 
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EXPERIMENTS ON CURRENT-DRIVEN THREE-DIMENSIONAL ION SOUND TURBULENCE 
PART I. RETURN-CURRENT LIMITED ELECTRON BEAM INJECTION 

R.L. Stenzel 

TRW Systems* Redondo Beach, CA 90278 
ABSTRACT 

The injection of an electron beam into a weakly collisions!, 
magnetized background plasma is investigated experimentally. The injected 
beam is energetic (40 eV) and cold (< 0.5 eV), the background plasma 
(n e a 10 9 cnf 3 , u p > <u c ) is initially isothermal (kT g = kT- = 0.2 eV). 

Beam and plasma dimensions are so large 43 cm diam, 250 cm length, 

Debye length <0.02 cm) that the system can be considered unbounded. 

The temporal and spatial evolution of the beam-plasma system is dominated 
by collective effects. High-frequency electrostatic instabilities 
(w = Up) rapidly thermal ize the beam and heat the background electrons 
(kT e » kT^)* The injected beam current is balanced by a return current 
consisting of background electrons drifting toward the beam source. The 
drift between electrons and ions [(kT e /m- )^ 2 « v d < (kT e /m e )^ 2 ] gives 
rise to an ion acoustic instability which develops into strong three- 
dimensional turbulence to be discussed in a companion paper. It is shown 
that the injected beam current is limited by the return current which is 
approximately given by the electron saturation current. Non-Maxwell i an 
electron distribution functions are observed. 


29 69 7-600H-RU-00 
Page B-3 


I, INTRODUCTION 

The injection of an electron beam into a stationary background plasma 
is a problem of general interest. Such situations arise, for example, in 
space plasmas (aurora), in laboratory plasmas (dc discharges) and in fusion 
type plasmas (relativistic E-beams). While collision dominated systems 
have been analyzed for many years in glow discharges'* collision! ess beam- 
plasma interactions are presently of main interest. Particular interest 
has been given to high frequency electrostatic beam-plasma instabilities 
near the electron plasma frequency t» which are effective in randomizing 

H 

2 3 

an initially cold, energetic electron beam. 5 The current closure in 
such a system is an important problem which has been given less attention 

4 

in small laboratory devices than in unbounded media. However, with the 
development of larger laboratory plasmas the return current problem has 
to be considered as an integral part of beam injection experiments. In 
the present investigation, it will be shown that the return current gives 
rise to strong low frequency electrostatic instabilities at frequencies 
up to the ion plasma frequency to - . Ion sound turbulence is effective in 

r * 

scattering background electrons which carry the current. Thus, the 
instabilities provide resistance to both the injected beam current and the 
return current. 

The present paper describes the properties of the electron beam and 
background plasma whereas the instability analysis and wave dynamics are 

C 

given in a companion paper . The experimental arrangement is kept as 
simple as possible so that the current flow can be readily accounted for. 

An electron beam is extracted from a uniformly emissive, large diameter (43 cm) 
cathode and closely spaced ( y 1 cm) highly transparent (71%) anode. The beam 
is injected along a uniform magnetic field (B 0 = 75G) into a quiescent after- 
glow plasma of comparable diameter produced by the beam source previously. 
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The plasma is not in contact with the chamber walls except at the floating 
end section. Net current flows only between the closely spaced anode and 
cathode. The beam current injected through the gridded anode is balanced 
by a return current drawn from the background plasma. The incident beam 
is not reflected but scattered in velocity space even before reaching the 
end of the device. It is shown that the return current and thereby the in- 
j’ected beam current are essentially limited to the electron saturation 
current drawn by the anode from the heated background plasma. The strong 
currant flow changes the electron distribution function. Skewed, shifted 
electron distributions are observed. 

The first part of the paper gives a description of the experimental 

setup and techniques. In Section II, the beam and plasma properties are 
presented and discussed. A brief summary of the observations is given in 
Section III since further results are presented in Ref. 6. 

II. EXPERIMENTAL ARRANGEMENT AND MEASUREMENT TECHNIQUES 

The experiment is performed in a large diameter, uniform, weakly 
collisional magnetoplasma. The device (Fig. 1) which has been described 
previously in connection with whistler wave experiments^ ^consists briefly 
of a low pressure (2 x 1(T 4 Torr krypton) dc discharge between a solid, 
oxide-coated cathode of 43 cm in diameter and closely spaced 1cm), 
highly transparent (71%) mesh anode. The plasma is confined by a uniform 
axial magnetic field (B Q - 75G) and terminated after 2,5 m by a floating 
permanent magnet end wall]*"* The plasma is not in contact with the side 
chamber walls. As shown in Fig, 2, the discharge is pulsed repetitively 
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(t = 25 msec). The first discharge pulse produces an afterglow back- 
ground plasma (10 9 < n Q < 10^ cm - ^, kT = kT- - 0.25 eV). At a desired 

“SGI 

afterglow time (t - 20 msec) a second independent beam voltage pulse is 

d 

applied to the cathode and the electron beam is formed by accelerating 
electrons in the cathode-plasma sheath of several Debye lengths in thickness 
(< 2 mm). The highly transparent anode is well outside the sheath. The 
plasma potential is close to that of the grounded anode since the dominant 
electron loss is at the anode. All measurements are sampled so that one 
can observe the temporal evolution of the beam-plasma system. Spatial 
information is obtained from a variety of movable probes shown in Fig. 1. 

Plane, one-sided Langmuir probes are located at the anode (No. 1), at 
z : 7 cm from the cathode (No. 2), in the middle of the device (No. 3) 

'i 

♦“i 

and on an axially movable shaft (No. 4). Two pairs of small cylindrical > 

Langmuir probes (0.1 mm diam x 2.5 mm length) movable in orthogonal 

0 

planes are used to detect and/or excite ion acoustic waves. An array probe 
is employed for multichannel simultaneous fluctuation measurements. There 
are additional high-frequency probes (not shown in Fig. 1) which are 
used to excite whistler waves, electron Bernstein and Landau waves for 
diagnostic purpose. 

III. EXPERIMENTAL RE5ULTS 

First, the background plasma properties are analyzed. Fig. 3 shows the 
density vs. afterglow time. At high densities (tu » w ) the absolute value 

pc 

8 11 

of n is reliably found from whistler wave dispersion measurements. 1 

12 13 

In the range > ud c resonance cones are useful for density diagnostics. * 

Although cylindrical Langmuir probes are reliable for relative density 

14 

measurements an absolute density determination from orbital motion theory 
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leads to uncertainties since the ion temperature is not accurately known. 

The ion saturation current is thus "calibrated" by absolute density values 
derived from wave dispersion measurements for electron plasma waves 
[Fig. 4(a)] and electron cyclotron harmonic waves [Fig. 4(b)]. The 
observed density decay is due to ambipolar diffusion along and 
surface recombination at anode, cathode and floating end end wall. 

The electron temperature decays from kT g = 2 eV during the discharge to 
kT g = 0.2 eV in the late afterglow (t fl > 15 msec). From the damping of 
ion acoustic waves the ion temperature is estimated to be kT^ < 0.2 eV. 

Further parameters of the background plasma are summarized in Table 1. 

The beam v^ cage pulse is appV _ . at an afterglow time t fl which 
determines the desired background plasma density. After injection beam 
and background plasma strongly interact and change their initial properties 
to new steady-state values. Fig. 5 shows Langmuir probe (No. 2) traces 
prior to beam injection (top) and at t $ = 10 psec after turn-on of the 
beam pulse (bottom). Prior to beam injection, there is essentially no 
difference in the traces with the probe normal pointing toward the cathode 
( <j) = 0) or in the opposite direction. After injection, the energetic 
beam electrons (V fa = 46V) are observed for <J> = 0 but not for 4> = 180°, 
indicating that there are no reflected beam electrons. The background 
electron flux (total current minus beam current for $ = 0) is about twice 
as large toward the anode (4> = 180°) than in the direction of beam injection 
where beam and background saturation currents are approximately equal . 

Thus the background plasma carries a return current equal to the injected 
current. The total saturation currents for <f> = 0 and 4>= 180° are equal, 
i.e. the plasma is free of a net current, consistent with the fact that the 
applied current flows only within the narrow cathode-anode gap. 
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A closer inspection of the traces reveals first that the saturation 
current ("knee") for <j> = 180° occurs at a lower potential than for = 0° 
which is characteristic for a drift in the background plasma. Second, 
the electron "temperature" for $ = 180° is about twice as large as for 
$ “ 0°, which indicates drift and possible asymmetry in the distribution 
function. By differentiating the plane probe characteristics and 
displaying dI p /dV p vs. velocity, the parallel electron distribution 
function is obtained (Fig. 6). It shows the cold fast beam and the skewed, 
slightly shifted background distribution. 

The drift velocity measurement from Langmuir probe traces is unfortunately 
not very accurate. One can only estimate that the drift velocity exceeds 
sound speed (c g = 2.3 x 10^cm/sec) by one to two orders of magnitude yet 
the drift is about one order of magnitude smaller than the electron thermal 

g 

velocity (v = 1.4 x 10 cm/sec). Conventionally, the drift velocity is 
determined from the electron flux and density (j = nev d ). 7he density is 
well known (n = 2 x 10 cm" at t = 24msec), the electron flux (j = 3mA/ 

w a 

a?) i s calculated from the measured cathode current (6A), surface area 
(1450 cm ) and anode transparency (71 %) . The calculated drift velocity 
is v d = 0.94 x 10^ cm/s ec {- 40 c s = .08 v 0 ). However, only if the 
distribution remains Maxwellian the drift corresponds to an actual shift 
in the distribution function. In the present case, the current is partly 
carried by the anisotropy in the distribution and partly by an actual 
displacement of the distribution function maximum relative to the essentially 
stationary ion distribution. Only the shift is relevant for the ion 
acoustic instability. Since v d »c g the threshold for the instability is 
well exceeded. The anisotropy could arise from the fact that the anode 
which draws electron saturation current (see Fig. 9) represents an 
absorbing boundary for electrons. Electrons heated by beam-plasma 
interactions in the interior of the plasma stream toward the anode and 
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only a fraction 50%) is reflected by the anode-cathode boundary. If 

the anode transparency were reduced, the distribution would approach 

1 5 

a half-Maxwellian. The classical scale length for Maxwellization 
(4 = v e /v ge * 3 x 10 3 cm) is large compared to the present distance of 
observation from the anode. 

The spatial properties of beam and background plasma have been 
investigated. Across the magnetic field beam and plasma parameters are 
uniform to within 10% over a distance of approx. 35 cm. 8,17 Along B^, the 
beam-plasma system changes on a scale length short compared to the co- 
llision mean free path. Fig. 7 shows Langmuir probe (No. 4; $ = 0) 
traces for different axial positions z from the cathode sampled at a 
fixed time t after beam injection. With increasing distance of propagation 
the beam broadens and the background electron temperature and density 
increase. This rapid scattering of the beam is well known to be caused 

*[ C 

by electrostatic instabilities. Although in an unbounded system a wide 
range of oblique whistlers 17 and Bernstein waves are unstable, the mode 
with' the shortest growth length is that near the electron plasma frequency. 
For a cold beam plasma interaction the predicted maximum spatial growth 
rate is gi^en by" 18 




,V3 



O) 


where and v b /v e are the ratio of beam to background plasma density 

and directed beam velocity to background electron thermal velocity, res- 

_2 

pectively. With the experimental conditions (n^A^ = 5.8 x 10 , 
v b / v e = 3.2, k r = 6,3 cm" 1 ) the theoretical growth length is k^" 1 = 0.4 cm. 
The beam distribution is observed to change after approx. 10 wavelengths or 
25 e- folding lengths from the injection point. The beam broadens, forms a 
plateau and finally, for z a 20 cm, forms a tail of a distribution with 
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continously decreasing slope, - 3f e /3 Vj j < 0 . This stable distribution is 
evident in Fig. 8 which shows probe (No. 3} traces in the middle of the 
device (z = 125 cm). Although the fast particle flux is still different 
for <f> = 0 and 4 = 180°, its magnitude is small compared to the background 
electron saturation current. The drift in the background distribution is 
net observable within measurement accuracy. The beam electrons have 
transferred most of their energy to plasma waves which, via Landau damping 
heat the background electrons. Since the anode is a strong energy sink, a 
steady-state temperature is established. At the far end of the device 
(z = 250 cm), the asymmetry in the energetic electron flux amounts to 
only a few percent of the total electron flux. Reflections of energetic 
electrons thus play an insignificant role in the current balance. Corres- 
ponding to the axial gradient of the injected beam current there is a 
gradient in the return current. When the end plate instead of the grid 
is used as an anode, the discharge current is also limited to the back- 
ground electron saturation current. However, this mode of operation is 
avoided because it exhibits two undesirable features. First, the discharge 
cannot be initiated at normal discharge voltages (V^ = 40 V). Second, 
after the discharge has been triggered, the current exhibits large amplitude 
( 6 I/I = 75%) low frequency (f * 1 kHz) fluctuations. The peak anode current 
is comparable to the value for the gridded anode which is free of fluc- 
tuations. 

The temporal behavior of the beam-plasma system has been investigated 
next. Figure 9 shows Langmuir probe (No. 1) traces at different times 
before, during and after beam injection. The probe traces exhibit the 
background plasma properties at the anode {4 = 180°, z = 0). Prior to 
beam injection, the plasma potential has a small positive value consistent 
with the low temperature of the afterglow plasma. During beam injection. 
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the knee in the probe trace occurs at a negative potential with respect to 
the grounded anode. Thus, the anode draws a current in excess of the 
electron saturation current from the background plasma. A comparison 
between Figure 9 and Figures 5 or 7 indicates that the knee in the I-V 
characteristics shifts from V p = -1 .5V at z - 0 to V p = + 2V at z = 7 cm. 

This shift indicates that the background plasma drift increases toward 
the anode. It does not mean that the plasma potential decreases toward 
the anode which would imply an electric field {<a V p /Az > = 0.5 V/ cm) 
opposing' the current. Langmuir probe data are not sufficiently accurate 
to measure the small dc electric fields in a collision! ess plasma. In 
the present case, the theoretical electric field based on classical re- 
sistivity and a current equal to the electron saturation current is given 

by E = v_ n v„/( 2/? e/m Q ) = 8 mV/cm. Although in the presence of strong 

19 

ion sound turbulence the resistivity is expected to be much larger, it 
appears unlikely to be directly measurable with Langmuir probes, especially 
in the presence of drift velocity gradients. Estimates of the electric 
field from variations in the local floating potential are subject to 
additional errors since the tail in the distribution function may vary 
with position. 

After termination of the beam injection pulse at t = 200 psec, the 
return current also vanished. The plasma potential becomes positive 
(7^= + 5 kT e /e) such that the electron current to the anode is reduced to 
the ion saturation current. Since high energy electrons are lost pre- 
ferentially the electron temperature and plasma potential decrease in 
time. 

The temporal behavior -of the beam current and voltage are shown in 
Figure 10. Two time scales can be distinguished, (i) For t s 40 psec 
(top traces) density changes due to ionization are negligibly small. 
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The beam current is essentially constant except for a turn-on transient 
pulse. The anode current is the sum of the background electron saturation 
current, I = A.n_e{kT Q /2^m_}^^, and a fraction (1-T) of the total beam 
or cathode current (T is the anode transparency). Since anode and cathode 
current are equal their value is given by I d = I /T. For the present 

g _ o n 

parameters (n = 7 x 10 cm , kT - 3 eV, anode area A = - 421 cm' } the 
e e a 

expected current is - 19A which agrees with the observed value 1^ 20A 

from Figure 10. (ii) For t > 40 ysec the density increases due to 
ionization. This allows the beam current to rise even though the beam 
voltage decays which is the result of discharging the supply capacitor 
(C = 2 mF). 

The beam current increases due to both ionization and heating. The 
temporal behavior of electron density and temperature has been measured 
in two cases. Figure 11(a) shows the density increase after turn-on of a 
40 V discharge voltage pulse. A high background plasma density is chosen 
so that reliable density data from propagation of whistler waves can be 
obtained. The initial density increase is exponential with a time constant 
t = 150 usee. With an ionization rate v- = 1.2 x 10 6 sec” 1 (Table 1) 
the relative concentration of ionizing beam electrons in the center of the 
device is n./n - (tv.) 1 - 0.6% which is about an order of magnitude 
smaller than in the beam injection region. The beam electrons lose their 
energy much faster by collective effects than by collision phenomena. 

After the initial rise the density approaches a steady-state value, i.e. 
particle production and loss by diffusion balance. During beam injection, 
the loss rate is higher than in the. low temperature quiescent afterglow. 

Figure 11(b) shows the electron temperature dependence during beam injection 
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at a lower background density (n g = 2 x 10^ cnT^). The temperature is 
derived from the measured sound speed of test ion acoustic waves. During 
the first 50 psec, the test wave damping is very strong hence the measurement 
accuracy is low. However, from Langmuir probe data it is known that there 
is a fast 10 psec) initial temperature rise (&kT fi >1 eV). The initial 
heating process is due to resonant wave-particle interactions with the elec- 
trostatic instabilities since collisional interactions between beam and 
background electrons occur on a slower time scale (see Table 1). Collisional 
heating may explain the subsequent slow temperature rise which in turn 
is limited by heat losses to the anode. 

Ionization phenomena can, of course, be avoided when the beam energy 
is kept below the ionization energy of the background gas (14 eV). Figure 
12 shows Langmuir probe (No. 2, $ = 0°) traces at different beam voltages 
both below and above ionization potential. For beam voltages V b < 20V 
[Fig. 12(a)] beam and background flux are almost independent of beam energy 
whereas for higher voltages [Fig, 12(b)] the currents rise with beam energy. 

The beam flux is always linked to the background flux and it is not possible 
to raise the relative beam current by increasing the beam voltage. 

The current- voltage characteristics of the beam-plasma system are 
summarized in Fig. 13. The characteristics depend on time since, in general, 
background temperature and density are affected by the beam. Steady-state 
conditions are approximately achieved at t s 0.5 msec. The characteristics 
are, of course, functions of initial density, magnetic field, gas, pressure 
and geometry, too. the present beam currents are far below the emission 
capabilities of the cathode. At higher pressures discharge currents of 
500A are readily obtained, and on smaller similar cathodes in strong 
magnetic fields current densities of j^=l A/cm are commonly observed. 
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IV. SUMMARY AND CONCLUSIONS 

Time and space resolved measurements of the basic parameters in a pulsed 

electron beam injection experiment have been presented. It is shown that 

in the essentially unbounded, uniform beam-plasma system the current of the 

injected beam closes via a field-aligned return current of background 

electrons. The maximum beam current which can be injected does not 

significantly exceed the background electron saturation current. High 

frequency electrostatic instabilities inhibit the penetration of the electron 

beam; the beam is scattered in velocity space on a distance short compared 

to the collisional mean free path. Low frequency instabilities limit the 

drift of the background plasma to a value small compared to the electron 

thermal speed. The analysis of the' high frequency instabilities has been 
17 

presented earlier, the low frequency instabilities will be discussed 
in a subsequent paper. 6 
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TABLE 1: BASIC PLASMA PARAMETERS 


Gas, pressure 

Krypton, 2.5 x 10“ 4 Torr 

Density (t a = 24 msec) 

n e : 

2 x IQ 9 cm" 3 

Electron Temperature 

kT e 

= 2 eV 

Ion Temperature 

kT i 

= 0.2 eV 

Magnetic Field 

V 

= 75 G 

Plasma frequency 

f pe 

= 402 MHz 


f pi 

= 1.05MHz 

Cyclotron frequency 


= 210 MHz 


f c1 

= 1.43 KHz 

Debye length 


= 0.23 mm 

Cyclotron radius 

r ce 

= 0.64 mm 


r ci 

=7.7 cm 

Electron collision frequency with 



neutrals 20 


s 4.2 x 10® sec" 1 

15 

ions 

electrons 1 ® 

en 

v ei 

v ee 

- 81 sec" 1 

= 3.1 x 10 4 sec" 1 

Electron mean free path 

*e 

= 200 cm 

Ion-neutral charge exchange 



collision frequency 20 

v in 

= 5 x 10 3 sec" 1 

Ion mean free path 

Beam electron mean free path 

(V. = 40V) for 
o 20 

elastic collisions 

£ i 

= 13 cm 


= 90 cm 

. . . , . . . 20 
ionizing collisions 

e 

*e 

= 320 cm 

20 

Ionization energy 

E i 

= 14 eV 
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FIGURE CAPTIONS 

1. Schematic diagram of the plasma device and diagnostic tools. 

2. Typical time sequence of the repetitively pulsed beam injection 
experiment (t^= 25 msec) the first discharge pulse creates an 
afterglow plasma into which the beam is injected at t . Sam- 
pling measurements show the temporal evolution of the beam-plasma 
interaction. 

3. Density decay of the afterglow plasma measured by independent 
techniques, B Q = 75 Gauss. 

4. (a) Dispersion measurements of test electron plasma waves propagating 
along if yield the density in the presence of a beam. 

(b) Dispersion measurements of electron cyclotron harmonic waves 
in the absence of beams. are used for independent density diagnostics. 

5. Langmuir probe traces in the absence (top 2 traces) and presence 
(bottom) of the injected beam. The one-sided plane probe faces 
the beam source at $ = 0° (see insert sketch) and the background, 
plasma at <f> - 180°. z = 7 cm, = 46 V, B Q = 75 G, n £ = 2 x lO^cnf^, 

6. Parallel electron distribution function derived by differentiating 
the bottom traces in Fig, 5. 

7. Plane Langmuir probe traces sampled at a fixed time at different 
distances z from the cathode. The traces show that beam and 
background heating occur over distances short compared to the 
electron mean free path { sl q - 90 cm). t fl = 24 msec, B q - 75G. 

8. Langmuir probe traces in the middle of the plasma device (z =125 cm) 
for the probe facing the beam ($ = 0) and opposing it (<j> = 180°). 

As a result of electrostatic instabilities the beam has been 
broadened and lost most of its energy on a scale length short 
compared to the ionization mean free path - 300 cm). 

c 
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Fig. 9. 


Fig. 10. 


Fig. 11. 


Fig. 12. 


Fig. 13. 


Plane Langmuir probe traces sampled at different times before, 
during and after the beam injection pulse. The probe is located 
at the anode (z = 0) and faces the background plasma ($= 180°). 
During beam injection the grounded anode draws currents in excess 
of the electron saturation current; in the absence of the beam 
the anode is at the floating potential. = 45V, t a = 24 msec, 

B q = 75 G, Kr 2 x 10“^ Torr. 

Cathode current and voltage vs. time on different scales ($ = common 
zero line). At early times (t < 40 nsec) the beam current is 
constant. It is limited by the return current carried by the back- 
ground plasma whose density increases on the longer time scale due 
to ionization. 

{a} Density increase during injection of an ionizing electron beam 
into a dense partially ionized background plasma. The density is 
determined from the dispersion of whistler waves. 

(b) Electron temperature vs. time during beam injection into a low 
density background plasma. The temperature is obtained from ion 
sound speed measurements; the dashed initial rapid rise is inferred 
from probe measurements. 

Langmuir probe traces at a fixed sampling time t and position 
(z = 7 cm, 0 = 0 °) for different beam voltages below {a) and above (b) 
ionization potential (14 eV). The beam current does not increase 
with voltage unless the background electron density and temperature 
are raised. 

Current- voltage characteristics of the beam source. Since the 

injected beam current is limited by the return current the 

characteristics depend on initial background plasma parameters and 

subsequent temperature and density changes due to the beam-plasma 
interaction. 
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FIGURE 7 
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FIGURE 11 
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EXPERIMENTS ON CURRENT-DRIVEN THREE-DIMENSIONAL ION SOUND TURBULENCE 

PART II. WAVE DYNAMICS 
R. L. Stenzel* 

TRW Systems, Redondo Beach, California 90278 
ABSTRACT 

A fast (40 eV), cold (< 1 eV), large diameter (43 cm) electron beam 

is injected along B Q (75 6) into a uniform background plasma (10 9 ■< n e < 

10 -3 

10 cm , kT « 10 kT. = 2 eV). The injected beam current is balanced 

C 1 

by a field-aligned return current which gives rise to an ion acoustic 
instability {c g v d < v , where c s , v d , v £ are sound, drift and electron 
thermal velocity, resp.). The plasma is essentially unbounded (charac- 
teristic dimensions 10 Debye lengths) and collisionless on the time scale 
of the pulsed experiment {v^ n /co ^ = 5 x 10~^). The fluctuations are 
analyzed by probe techniques in real time and space as well as in the 
frequency and wavenumber domains. Strong density fluctuations are ob- 
served in the entire ion acoustic spectrum reaching peak amplitudes 
( 5n e / n e ~ 20%) at the low frequency end (w/w . =0.2). The shape of the 
frequency spectrum is close to that predicted by renormalized plasma 
turbulence theory. The current is pulsed and the development of the 
instability from linear growth to saturation is observed. Spatial 
correlation measurements indicate that the turbulent wave vector spectrum 
extends in magnitude to kx^ > 1 and in angle up to ± 90° with respect 
to the electron drift direction. The propagation of phase coherent test 


* Also at Department of Physics, University of California, Los Angeles, 


California 90024. 
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waves in the current carrying plasma is investigated. For the highly 
turbulent regime strong damping is observed while at lower drift velo- 
cities well defined wave propagation within a half space centered around 
the drift velocity direction is observed. Phase and amplitude contours 
of test waves in orthogonal planes are presented. Oblique propagation 
of two strong test waves at different frequencies indicate the role 
of nonlinear ion Landau damping in absorption of wave energy. 
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I . INTRODUCTION 

The current-driven ion acoustic instability is a well known pheno- 

1-6 7 8 

tnenon which has been analyzed in theory , computer simulations 5 , and 
9-1 ft 

many experiments. While the linear properties are understood the non- 
linear saturation which may involve strong turbulence remains to be ex- 
plained. Experimentally, the instability is produced either by driving 
a toroidal current with an induced electric field ^ or by drawing a 
current between two electrodes. The latter approach has been thoroughly 
studied in positive col umns^ 5 ^’^’^, i.e. collisional, bounded and 
unmagnetized plasmas. Some work has been done in weakly collisional 
plasmas' ,0 ’ 12,13,15 ’' l6 few of which can be considered free of boundary 
effects . 

In the present experiment a current driven ion sound instability is 

observed in an essentially coilisionless and unbounded plasma immersed in 

a uniform magnetic field. The current corresponds to the field-aligned 

return current associated with the injection of a fast electron beam into 

19 

a magnetized background plasma as described in the companion paper. 

The current is pulsed and the evolution of the instability is observed. 

The instability saturates within approximately 15 ion plasma periods at 
a level of strong turbulence (W/n kT = 5%, where W, n kT are wave energy 
and kinetic particle energy, respectively). Distribution function 
measurements^ show that the drift velocity is much larger than sound speed 
although it remains well below the electron thermal velocity (v^ < 0.1 v g ). 

The distribution function is not Maxwellian. Analysis of the unstable 
wives indicates that the propagation is 4 in general, three dimensional. 

The correlation lengths are approximately 50 x d (x d = Debye length) in the 
highly turbulent state. Simultaneous observations of density fluctuations 
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in time at different locations with an array probe show the microscopic 
properties of turbulence, i.e. growth, propagation and decay of pertur- 
bations . 

The propagation of test ion acoustic waves in the current carrying 

plasma gives additional insight into the instability properties- While 

1 5 

previous investigations were restricted to one-dimensional propagation 
along the electron drift the present experiments show the three-dimensional 
propagation characteristics. At lower turbulence levels weakly damped test 
wave propagation is observed in an angular half space centered about the 
direction of the drift velocity. Propagation of two oblique test waves of 
large amplitudes results in an off-resonant beat mode which is strongly 
Landau damped when the phase velocity is near the ion thermal speed. Non- 
linear ion Landau damping has been proposed as an important mechanism for 
energy absorption and transfer among different modes. 

Since the experimental setup and techniques have been described 
19 

earlier the present paper starts with the measurement results. The 
characteristics of the instability are given in Section II and the test wave 
experiments are presented in Section III. The results are discussed and 
summarized in the final Section IV. 

II. INSTABILITY ANALYSIS 

The presence of a low frequency instability is readily inferred 
from the large fluctuation level in the Langmuir probe saturation current. 

While the previously shown probe traces^ represented time-averaged values 
Fig. 1 displays a fast sweep of the probe characteristics during the pulsed 
beam injection. The probe voltage (top trace) varies linearly in time; the 
probe current (bottom trace) is recorded for three single shot exposures. 

One can see large amplitude (sn_/n Q = 25%), low frequency (f - 200 kHz) 

6 S 
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fluctuations in the electron saturation current which are uncorrelated 
from pulse to pulse. The current fluctuations are not due to fluctuations 
in the plasma potential but represent density fluctuations. When a high 
impedance probe circuit is used the fluctuations are also observed in the 
ion saturation current. However, the frequency response is found to be 
much better in the electron saturation regime than for iori collection. 

The presence of strong density fluctuations is confirmed by probe 

independent diagnostics. A resonance cone is excited with a short-wire 

antenna at a frequency u < w c < Up. 1 Since the cone angle is a function 

of density low frequency density fluctuations scatter the cone angle. When 

the rf field is detected with a second antenna the received signal is 

amplitude and phase modulated as a result of the spatial deflection of 

22 

the cone with respect to the fixed antenna. Fig. 2 shows the frequency 
spectrum of a small amplitude resonance cone signal in the presence of the 
return current. The lower sideband indicates the presence of broadband 
density fluctuations. The shape of the modulation spectrum varies with 
probe location across the cone. The sideband vanishes in the absence of 
the return current. Although this observation is rather qualitative it 
nevertheless confirms that the density perturbations detected with Lang- 
muir probes are real and not caused by probe perturbations on the plasma. 

Due to the pulsed nature of the experiment growth and saturation of 
the instability can be observed. Fig. 3 shows on the top trace the pulsed 
beam voltage {V^ « 35 V) vs time and on the bottom trace the electron 
saturation current to a small cylindrical Langmuir probe (0.13 mm diam., 

2,5 mm length). Several traces are superimposed so as to indicate the 
transition from coherence to turbulence. The instability grows coherently 
for approximately two periods of oscillations at f * 400 kHz and then goes 
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over into turbulence which is fully developed and saturated (6l e /I e = 
constant) after 10 to 15 ion plasma periods (f * = 1 MHz). The increase 

r * 

in the average electron saturation current is due to heating of background 
and beam electrons and the associated change in the saturation charac- 
teristics of small probes! The observed growth characteristics are 
essentially independent of probe position* thus the instability starts 
throughout the plasma volume. Since the coherent growth time is short 
the subsequent observations will apply to the saturated state. 

The frequency spectrum of the saturated instability is shown in Fig. 4. 

The noise extends through the entire ion acoustic frequency range but 
exhibits a pronounced maximum { 6n / n = 20%) at the low frequency end 
(tu/tOp^ = 0.2). The low frequency peak in the spectrum is not the result 
of collisions (v. /w . = 0.005), ion cyclotron modes (u -/w • = 0.001) or 

i n pi ^ i pi 

boundary effects (c s /D = 5 kHz, D is the plasma diameter) but is a charac- 
teristic property of the stationary turbulence spectrum. According to 
renormalized turbulence theory which considers both scattering of ions and 
electrons by the waves the long wavelength Landau resonance is destroyed by 
electron scattering while the short wavelength intensity fall-off is the 
result of nonlinear ion Landau damping. The observed spectrum has the 
same characteristic shape as that predicted by Horton et al . (See Figure 2 in 
Ref. 6). The spectrum may be considered stationary since its shape does 
not significantly vary in time provided the slowly rising value of is 
taken into account. The wave energy normalized to the particle energy 
density is W/n fi kT e ** (sn e /n 0 ) 2 = 5x 10" 2 . This is a large value; for 
example, the wave electric field at the peak noise frequency f = 200 kHz 
(X = 5 mm) is 1 6E j = k5p = 2.5 V/cm. Although 6$ decreases with increasing 
to the electric field remains large as to ■+ w 
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The wavenumber spectrum is determined by correlation techniques. The 
broadband signals (frequency range 0.05 - 5 MHz) of two identical, movable, 
cylindrical Langmuir probes (0.13 mm diam, 2.5 mm length, 1 mm diam shield) 
are multiplied, sampled and time-averaged. Either probe signal ^ , $ 2 can be 
delayed by a broadband delay line. For three-dimensional wave propagation 
the measured unnormalized cross correlation function (Ar, x) = J'i- |(r,t ) 
$ 2 (r + Ar, t + t) dt is a function of three spatial coordinates Ar = (ax. Ay, 
az) and delay time t. First, in Fig, 5(a), the correlation signal is 
displayed vs one variable, the coordinate Az along the electron drift. The 
half width of the correlation signal, az * 1 mm, gives a measure for the 
spread in the magnitude of the wavenumber k j j . It does, however, not indicate 
the spread in the angular k-vector spectrum. Fig. 5(b) displays the cross 
correlation function vs two coordinates along and across B Q . The value of 
C^ 2 (0> Ay, az, 0) is indicated by contours of constant relative amplitude. 

The observed elliptical shape of the correlation surface indicates a 
strong angular spread in wave propagation with respect to the drift direction 
along B 0> 

The direction and velocity of wave propagation is determined by dis- 
playing C^ 2 as a function of delay time t and position variables. Fig. 6 
shows the propagation characteristics along B 0 . With increasing delay 
(t > 0) of the movable probe signal (No. 2) the peak of the correlation 
signal shifts to negative values of Az, and vice versa for delay (t < 0) 
in the stationary probe signal (No. 1, located at Az = 0). Thus, the wave 
propagation points in the direction of the electron drift. From the 
diagram of Fig. 6(b) the parallel phase velocity is found to be v 0 = 

1,9 x 10 cm/sec which corresponds to sound speed in Krypton at kT g = 3 eV. 
Further useful information can be obtained from the spatial amplitude 
decay of the cross correlation signal. The e-folding decay length which 
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is a measure of the coherence length of the waves is found to be 
az c = 4 mm = 50x d . 

The two-dimensional cross correlation function shown in Fig. 5(b) 
does not completely describe the oblique propagation characteristics. 

The frequent assumption of propagation symmetry in the plane normal to 
the drift direction is not supported by the present observations. Fig. 

7 shows cross correlation measurements in a plane normal to B . The 
correlation signal along the coordinate Ax normal to the plane of Fig. 

5(b) is shown in the upper part of the figure. By taking a number of 
such traces at different positions Ay, except for Ay = 0 where probe 
shafts interfered, the map of Fig. 7(b) is obtained. Although no gradients 
in the bulk plasma properties have been noticed there are preferred direc- 
tions of wave propagation normal to the major axes of the elliptical con- 
tours. Whether this asymmetry is the localized wave propagation in £-space 

O 

proposed by Kadomtsev or is a result of minor gradients in the electron 
drift direction is difficult to assess. The minimum half widths of the 
correlation peaks are comparable to the values of the axial correlation 
function. Thus, in different planes (Ar, az) one can find varying degrees 
of angular spread around v^ ranging from 9 = arcos {v d /c g ) = 90° to nearly 
zero. 

The direction and speed of wave propagation in the perpendicular plane 

are indicated in Fig. 8. In contrast to Fig. 6(a), the correlation function 

does not shift with increasing delay in a unique direction but splits and 

* 

moves in opposing directions. Thus, there are waves propagating both to 

the right and to the left along the axis Ax, The propagation speed (v = 

5 p 
2.2 x 10 cm/sec) is larger than sound speed indicating that the waves are 

again oblique to the chosen coordinate. Clearly, for a complete description 
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of the wave properties the full cross correlation function depending 
on all 4 variables has to be measured. In principle, this can be done 
with the present techniques but, in practice, this formidable task 
requires the aid of a computer for data handling which was, unfortunately, 
not available. 

A different approach to characterize turbulence is the direct dis- 
play of the fluctuations in real time and space. This is done with an array 
probe and a multi-channel data display as shown in Fig. 9. The signals 
from six identical, radially displaced, cylindrical Langmuir probes are 
displayed simultaneously on three dual beam oscilloscopes. A single shot 
exposure of the fluctuations vs time at different points in space is shown 
in Fig. 10. The figure displays the microscopic properties of turbulence: 

There are density perturbations which grow (a) or decay (b) in the direction 
along the array. From the time delay and probe spacing one concludes that 

these are waves propagating close to sound speed (v = 1.4 x 10 5 ^-) in the 

P S6C 

direction from probe No. 6 to No. 1 essentially perpendicular to the elec- 
tron drift. The amplitude growth and decay thus reflect the coherence 
properties rather than oblique propagation effects. The typical "lifetime" 
of a density perturbation is 7 ysec during which it propagates about 
10 mm. There are also waves which propagate oblique to the array (c,d) as 
indicated by higher phase velocity components (v * 2.5 x 10 5 ~ ) and 

bet 

smaller amplitude variations following more along the crest of the waves. 

By expanding the array to two or three coordinates one can, in principle, 
determine the propagation direction, velocity and amplitude behavior of 
individual density fluctuations. 
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III. TEST WAVE PROPAGATION 

By launching small amplitude monochromatic test ion acoustic waves 

from a well defined source useful information about the growth rate and 

saturation process of the instability can be obtained. In the presence 

of strong turbulence the ion acoustic wave dispersion is expected to be 
23 24 

modified * , and test wave experiments on new ion sound modes have been 

pc pc ' 

reported. * In rather few experiments three-dimensional wave propa- 

27 

gation effects have been considered. 

In the present test wave studies the waves are excited near the anode 
(z = 7 cm) from one of the short cylindrical Langmuir probes (0.13 nun 
diam, 2.5 mm length) and received with the second identical probe. In the 
plane normal to the exciter wire azimuthal ly symmetric cylindrical waves 
are excited in the absence of drifts. The receiver probe can be scanned 
in this plane, and the plane can be inclined with respect to the magnetic 
field. Thus, with minimum perturbation on the plasma three-dimensional 
propagation effects can be investigated. 

Propagation of a Single Mode 

In the measurements presented here the exciter signal is a single 
frequency cw signal. The receiver signal is usually amplified and applied 
to a broadband interferometer circuit whose output is proportional to 
$>■] (r, t)cos (£(.{.). r) where ^ (r, t), £{r, t) are wave potential and - 
propagation vector, resp., which may vary in space and time slowly com- 
pared to a wavelength and oscillation period, respectively. Fig. 11(a) 
shows the interferometer signal vs time at a fixed frequency (f = 3 MHz) 
and axial probe spacing (az « 12 mm). Also indicated are the beam pulse 
(top) and the sampling pulse (bottom trace). The growing oscillations 
indicate that amplitude and wave number vary in time. The former is 
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mainly due to the density increase which raises to above to and 
then reduces ion Landau damping as to « u . , and partly due to the lower 

r 1 

level of turbulence at later times when the current penetration decreases. 
The latter is mainly due to electron heating {k = w/c s = T "^) and 
partly due to dispersion [k - k (tUp^ , T g ), see Eq. (1)]. The fluctuations 
in the traces indicate the scattering of the test wave by the background 
noise which is not directly visible due to some low pass filtering of the 
interferometer signal. At early times (t s < 100 ysec) test wave propa- 
gation at any frequency u < w . exhibits severe damping due to wave-wave 
interactions in the strong turbulence. After the end of the current 
pulse the amplitude sharply drops since the instability vanishes. The 
wavenumber increase due to electron cooling gives again rise to an os- 
cillatory pattern. 

The temporal variation of the wavenumber gives rise to a phase 
modulation of the test wave signal which is either observable in the time 
domain [Fig. 11(a)] or the frequency spectrum as shown in Fig. 12. In the 
latter case the received signal $^cos (wt - kAz) is directly applied to 
a spectrum analyzer whose output is sampled. During the beam pulse 
[Fig. 12(b)] the decrease in k produces an upshifted sideband, in the 
afterglow of the beam pulse [Fig. 12(c)] the increase in k lowers the fre- 
quency. If the frequency spectrum were the only data available as is 
common in many pulsed experiments the interpretation of the sidebands 
could be ambiguous since Doppler effects and wave scattering produce 
similar results. 

In Fig. 11(b) the interferometer signal is sampled at a fixed time 
t s and displayed vs probe separation Az. The test wave propagates at an 
almost constant amplitude in the direction of the electron drift. The 

C 

phase velocity (v^ = fXj| * 1.7 x 10 cm/sec) corresponds to sound speed 
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in Krypton at kT g = 3 eV* a value approximately found from Langmuir 
probes] 9 By varying the sampling time the change in wavelength or 
electron temperature is verified. 

At a fixed sampling time the test wave dispersion relation w vs 
k | | has been measured (Fig. 13). Due to the presence of the de- 
stabilizing electron drift waves almost up to the ion plasma frequency 
(f .. = 3 MHz) can be launched from a simple wire exciter. Since the 
probe dimensions are small (0.13 mm diam) and a high precision probe 
drive system has been used wavelengths as short as 300 um (kx^ = 1,3) 

have been measured which is commonly considered the regime restricted to 

1 8 

laser scattering diagnostics. The figure shows a constant phase velo- 
city equal to sound speed c g = X D at low frequencies and dispersion 
effects near the ion plasma frequency as predicted by the dispersion 
relation for ion acoustic waves from fluid theory 

kx D 

ti) = 01 • n 9 1/9 0 ) 

pl (i + k 2 x 2 ) 1/2 

The detailed amplitude and phase distribution of a monochromatic test 

wave in a plane along t Q .(coordinate z) and across B Q (coordinate y) 

* 

is shown in Fig. 14. The exciter probe is located at the origin with its 
axis normal to the measurement plane. Fig. 14(a) presents a raw-data 
interferometer trace along ^ which indicates weakly damped large ampli- 
tude ion waves propagating parallel to the electron drift v^ , and strongly 
decaying small amplitude waves in the direction opposing v^. Since the 
test wave is not a plane wave but diverges the constant amplitude pattern 
implies spatial growth. 


1 
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Fig. 14(b) shows the phase fronts, i.e., contours of the interfero- 
meter maxima (circles) and minima (crosses). There is a noticeable 
difference in the parallel wavelengths or phase velocities along 

r* _v C 

v d (v p+ = fx + = 1.86 x 10° cm/sec) and opposing v d (v p _ e 1.02 x 10 
cm/sec). A possible explanation for the asymmetry is a Doppler shift due 
to an ion drift toward the anode. The drift velocity is approximately 
v di = ^ v p+ " v p-^ 2 “ 4,2 x 1C>4 cm / sec > i * e * v di/ c 5 = 0-29 where c £ = 

(v . + v )/2 = 1.44 x 10 cm/sec. In a quiescent plasma the ions could 

p+ p- 

not possibly drift against the electric field associated with the elec- 
tron return current. However, in the present turbulent plasma where the 
wave electric field ( | sE | =2.5 V/cm) is much larger than the dc field 
(E dc = m v d v g /e <0.14 V/cm, where v d < 0.1 v g = 10 7 cm/sec and < 
u - = 2 tt x 4 MHz) the ion motion may be considerably different. The 
large amplitude waves can trap a significant number of ions resulting 
in asymmetric energetic tails or an average ion drift in the mean direction 
of wave propagation i.e,, parallel to the electron drift. Although 
direct ion distribution function measurements have not been performed in 
this experiment the ion tail formation has been observed in other experi- 

90 . 7 

ments and in two-dimensional computer simulations. 

Fig. 14(c) shows contours of constant relative wave amplitude. In 
a current-free plasma the expected contours are concentric circles around 
the exciter since phase and group velocities point into the same direction 
(uj >> to ,. ) . In the present case the amplitude contours are not concen- 
tric with the source but are displaced in the electron drift direction. 

Thus, for az > 0, the cross section ax * Ay of a constant amplitude region 
.increases initially with distance from the exciter which can be interpre- 
ted as wave amplification proportional to the increase in cross-section. A 
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quantitative evaluation of the plane wave growth rate would involve a 
calculation of the three-dimensional probe radiation pattern which has 
so far not been done. 

Under identical conditions the test wave propagation has been in- 
vestigated in the plane normal to the magnetic field and the results are 
shown in Fig, 15{a) - (c). Concentric circular phase fronts indicate 
the absence of ion drifts and electron temperature gradients in this 
plane. However, the amplitude contours are not completely symmetric and 
there are preferential directions for weak and strong amplitude decay. 

Similar asymmetries were noted in the noise correlation measurements and 
are possibly caused by minor gradients in the drift velocity. 

Propagation and Interaction of Two Modes 

While the experiments with a single ion acoustic test wave show basic 
dispersion and propagation characteristics the interaction among waves 
which will occur in the turbulent plasma can be investigated with multi- 
ple test waves of different frequencies. Three-dimensional wave-wave 
and wave-particle interactions are important mechanisms in establishing 
the stationary turbulence spectrum. ’ u In current-free, unmagnetized 

plasmas experiments on nonlinear ion Landau damping have been reported for 

♦ » 

one-dimensional ion acoustic wave propagation" and these have recently 
been extended to three dimensions. In the present work these effects 
are investigated in a current carrying plasma where the coupling coefficients 
are much larger than for v^ = 0. 

As shown in Fig, 16(a) two test waves of different frequencies are 
launched from separate cylindrical exciter probes at different radial 
positions. With a third identical Langmuir probe the phase and amplitude 
distribution of the two oblique test waves {u^, ^ , u 2 , ^ anc * the 


c-x* 
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nonlinear beat mode (w 3 = ^ - ugj - t-j - Fg) is mapped in the plane of 
the two exciter probes and \Tj. Interferometer traces parallel to v d along 
the symmetry axis of the two source are shown in Fig. 16(b). The ratio 
of the absolute wave amplitudes is $-j : $2 : $3 “ 1*2:1. Since the 
beat mode is destabilized by the electron drift its amplitude is comparable 
to that of the two test waves, i.e. about two orders of magnitude larger 
than in the current free case. However, since in the present case the 
beat frequency is near the peak of the turbulence spectrum whose saturation 
amplitude is comparable to the maximum beat mode amplitude the signal -to- 
noise ratio is of order unity. While the frequency matching condition for 
the beat mode is accurately satisfied the verification of the wave vector 
matching requires two-dimensional interferometer data. 

Fig. 17 displays the phase fronts of both excited modes cind the 
nonlinear beat mode. For oblique propagation the beat mode is a nonresonant, 
driven mode with phase velocity v p = U-| - o^)/)]^ - f 2 [ < U-j - w 2 )/ 

(|&j [ - |icg|) = c £ . The local wave vector matching condition is found to 
be well satisfied as indicated by the two examples shown in Fig. IS 
(points A and B are defined in Fig. 17). The wave vectors It-j , itg and ^ 
are determined from the local wavelength and wave normal direction. For 
comparison, the expected beat wave vector Alt = Itj - £g is shown together 
with the measured vector £ 3 and the two it values agree to within measure- 
ment accuracy. The difference in the direction of It,, and {£> changes the phase 
velocity of the beat mode w 3 /|lc 3 ! drastically. While at point A the beat 
mode still propagates near sound speed, at point B its phase velocity falls 
into the ion distribution furrtion and the wave is subject to strong ion 
Landau damping. The observed amplitude distribution is shown in Fig. 19. 

Contours of constant amplitudes analogous to Figs. 14(c), 15(c) have been 
measured for the two excited waves and, as displayed in Fig. 19(a), the 
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amplitude product ^ *$> 2 has been formed since this quantity is propor- 
tional to the driving term for the beat mode. The measured amplitude 
distribution of the beat mode is shown in Fig. 19(b). The beat mode does 
not peak near point B where ^ • $ 2 is largest because of its slow phase 
velocity in this region [v - (kT^/rn^ It rather maximizes near 

point A where ion Landau damping is small even though the driving term 
$1 • $ 2 is reduced. Although the importance of electron over ion non- 

Op C 

linear Landau damping has been pointed out * the present observation 

1 /2 

of strong beat mode absorption when v (kT./m-) indicates that ion 

r ■ * 

Landau damping cannot be neglected in three-dimensional large amplitude 
ion sound turbulence. 


IV. SUMMARY AND CONCLUSIONS 

An experimental investigation into low frequency instabilities 
associated with the pulsed injection of a fast electron beam into a 
magnetized background plasma has been presented. By analyzing frequency 
and wavenumber spectra as well as the electron distribution function the 
noise has been identified as ion acoustic waves driven unstable by the 
return current to the injected beam. The instability evolves rapidly 
into strong ion sound turbulence. The three-dimensional aspect of the 
wave propagation has been emphasized. The low frequency cut off of the 
instability predicted from renormalized plasma turbulence theory has been 
verified in this experiment which is essentially free from boundary 
and collision effects. The magnetic field effects of ions are negligible 
and no ion cyclotron modes or lower hybrid modes have been identified. 
Propagation of test waves confirm the basic features of the ion acoustic 
instability. Nonlinear wave-wave interactions between two large test 
waves lead to beat modes which in turn are unstable and build up to large 
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amplitudes. Multiple interactions among many such modes can give rise 
to the observed wave turbulence. If the beat modes interact with the 
ions strong wave damping is encountered. 

Besides the direct observations some conclusions may be drawn in- 
directly from the data. For example, the larger ion sound speed in the 
direction of the electron drift is indicative of an asymmetric ion dis- 
tribution created by ion trapping in large amplitude waves. Electron 
trapping may also occur since the bounce frequency = k{e^/m e )^ 2 = 

Q T 

4 x 10 sec {k * 12 cm , ^ - 0.5 V) is far above the highest electron 

C 1 

collision rate \> = 4 x 10 sec . However, since the electron Larmor 

en 

radius is comparable to the wavelength trapping would mainly occur parallel 
to B . Trapping may flatten the distribution function so as to saturate 
the instability , produce asymmetric distributions and reduce the current 
flow. Since highly oblique waves are present many electrons can resonate 
with these waves (v g =■ w/kj j where w < 0 ^., kj|/k^< ' ) but the 

nearly perpendicular wave electric field gives rise to only a small 
drift velocity (£/B Q * 3 x 10 6 cm/sec). 

Although numerous observations have been presented the analysis of the 

beam-plasma interactions is not yet complete. Direct measurements of the ion 

distribution and of dc electric fields still have to be performed. Further 

0*1 

more, the interaction between high frequency wayes^ driven unstable by the 
injected fast electron beam and the low frequency instabilities of the return 
current remain to be investigated. 

ACKNOWLEDGEMENTS 

The author wishes to thank Drs. R. W. Fredricks, W. Gekelman, W. 
l ! DiVergilio and M, Caponi for helpful discussions. The expert technical 

assistance from Mr. W, Daley is gratefully acknowledged. 



29697-6002-RU-00 
Page C-19 


This work was supported in part by TRW Independent Research and 
Development Funds and in part by NASA under Contract NASW-2953, 


29697-6002-RU-0Q 
Page C-20 


REFERENCES 

1. E. A. Jackson, Phys. Fluids 3_, 786 (1960) 

2. 8. D. Fried and R. W. Gould, Ph.ys. Fluids 4, 139 (1961 ) 

3. B. B. Kadomtsev, Plasma Turbulence (Academic, New York, 1965) p, 68 

4. L. M. Kovrishnykh, J. Expt. Theoret. Physics (USSR) 51_, 1975,(1966) 
[Soviet Phys. JETP 24, 1210 (1967)] 

5. J. A. Wesson and A. Sykes, Phys. Rev. Lett. 31_, 449 (1973) 

6. W. Horton, Jr., Duk-In Choi, and R. A. Koch, Phys. Rev. A. 1_4, 

424 (1976) 

7. D. Biskamp and R. Chodura, Phys. Rev. Lett. 27, 1553 (1971), Phys. 

Fluids 16, 893 (1973) 

8. M. Lampe, W. H. Manheimer, J. B. McBride, J. H. Orens, K. Papadopoulus, 

R. Shanny, and R. N. Sudan, Phys. Fluids 1_5, 652 (1972) 

9. V, Arunasalam and S. C. Brown, Phys. Rev. 140 A, 471 (1965) 

10. H. Tanaca, A. Hi rose, and M. Koganei, Phys. Rev. 161 , 94 (1967) 

11. S. M. Hamberger and J, Jancarik, Phys. Fluids 15., 825 (1972); s. Q. Mah, 
H. M. Skarsgard, and A. R. Strilchuk, Phys. Rev. Lett. 25_, 1409 (1970) 

12. Y. Amagishi , Phys. Rev. Lett. 29, 405 (1972) 

13. S. Watanabe, J. Phys. Soc. Japan 35, 600 (1973) 

14. D. B. Fenneman, M. Raether, and M. Yamada, Phys. Fluids 16., 871 (1973) 

15. A, Y. Wong, B. H, Quon, and B. H. Ripin, Phys. Rev. Lett. 1299 (1973) 

16. A. Mase and T. Tsukishima, Phys. Fluids 18., 464 (1975) 

D. B. II ic, Phys. Rev. Lett. 34, 464 (1975) 


17. 



29697-6002-RU-00 
Page C-21 

18. R. E. Slusher, C. H. Surko, D. R. Moler, and M. Porkolab, Phys. Rev. 

Lett. 36, 674 {1976) 

19. R. Stenzel, Phys. Fluids (companion paper) 

20. R. K. Fisher and R. W. Gould, Phys. Fluids 1_4, 857 (1971) 

21. R. Stenzel, Bull. Am. Phys. Soc. 21_, 1176 (1976); also submitted to 
J, Geophys. Res. 

22. R. Stenzel and W. Gekelman, Phys. Fluids, Jan. 1976 

23. S. Ichimaru, Phys. Rev. 165 , 251 (1968) 

24. T. Tanaka and K. Nishikawa, J. Phys. Soc. Japan 30, 894 (1971) 

* 

25. S. Watanabe and H, Tanaka, Phys. Lett. 33A, 125 (1970) 

26. Y. Kawai and M. Hasegawa, J. Phys. Soc. Japan 40., 607 (1976) 

27. T. Mantei , Bull. Am, Phys. Soc. 21_, 1076 (1976) 

28. R. D. Bengtson, K. W. Gentle, J. Jancarik, S. S. Medley, P. Nielson, 
and P. Phillips, Phys. Fluids 18, 710 (1975) 

29. T. Ohnuma, H. Izawa, and S. Adachi; Phys. Rev. A 12., 1648 (1975) 

30. V. N, Tsytovich, lisp. Fiz. Nauk 90., 435 (1966) [Sov. Phys. Uspekhi 
9, 805 (1967)] 

31. H. Ikezi and Y, Kiwamoto, Phys. Rev. Lett. 27, 718 (1971) 

32. M, L. Sloan and W. E. Drummond, Phys. Fluids 13., 2554 (1970) 

33. K, Nishikawa and C. S. Wu, Phys. Rev. Lett. 23., 1020 (1969) 


Z9697-6002-RU-00 
Page C-22 


Fig. 1. 


Fig. 2. 


Fig. 3. 


Fig. 4. 


o 


FIGURE CAPTIONS 

Fast sweeps of Langmuir probe voltage V (top trace, 5 V/div.) 

r 

and current I (bottom trace, 0.5 mA/div.) during the injection 

r 

of an electron beam (V^ = 10 V, t g = 100 vsec). Large amplitude, 

incoherent fluctuations in the electron saturation current are the 

result of the return current driven ion sound instability. The 

plane probe (No. 2 in Fig. 1, Ref. 19) is located at a distance 

az - 7 cm from the anode with its surface normal parallel to the 

-4 

injected beam velocity. Krypton 2 x 10 Torr, t = 24 msec, 

d 

B q = 75 G. 

Frequency spectrum of an rf resonance cone signal excited in the 
presence of the ion acoustic instability. Density fluctuations 
modulate the cone angle and generate the sideband to the applied 
monochromatic rf signal at w 0 / 2ir = 90 MHz. 

Top trace: beam voltage (20 V/div.); bottom trace: electron 

saturation current to a cylindrical Langmuir probe located at 
Az = 7 cm from the anode. Superposition of 8 exposures shows the 
development of the instability from coherent growth to turbulence 
at t = 15 ysec after turn-on of current. Probe bias +12 V, 

= 24 msec, Krypton 2.3 x 10"^ Torr. 

a 

Frequency spectrum of the saturated current-driven ion acoustic 
instability. The noise amplitude is the fluctuating electron 
saturation current of a cylindrical Langmuir probe Fourier trans- 
formed in a spectrum analyzer and sampled at t g . The absolute 
peak density fluctuation level at a/oip. = kA^ = 0.2 is 6n e / n e £ 

20%. The high frequency fall-off is fin Q /n_ * f~^. The noise 

3 

peak is not determined by collisions (v. = 5 x 10 sec ) or 
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finite geometry but scattering of electron by waves. Krypton 
2.3 x 10~^ Torr, z = 7 cm. 

Fig. 5(a) Unnormalized cross correlation function for delay time t = 0 
vs. position az along, B Q . The half width Uz^g = 1 mm) gives 
a measure for the width of the wavenumber spectrum, 

(b) Contours of constant correlation signal for t = 0 in the hori- 
zontal plane Az \ | B Q , Ay J_ B Q . The stationary reference probe ■ 
is located at Ay = az = 0. Elliptical contours indicate angular 
spread of wave propagation with respect to the preferred 
direction parallel to v^. Subsequent data (Fig. 7) show that 
propagation out of the plane also exhibits preferred directions, 
t = 24 msec, z = 7 cm, Kr, 2.3 x 10”^ Torr. 

a 

Fig. 6(a) Cross correlation signal at different delay times t vs. position 
az (ax. Ay = const.) The shift of the correlation peaks with 
delay time indicates a unique direction of wave propagation. 

The amplitude decay gives a measure for the correlation length 
of waves (az c « 5Qx D ) . 

(b) Time-position diagram for the correlation function (maximum A, 
minimum B), Sign and magnitude of the trajectory slope show that 
waves propagate at sound speed parallel to v^. 

Fig. 7(a) Cross correlation signal vs. position Ax across l£ normal to 
the horizontal plane of Fig. 5(b). 

(b) Contours of constant cross correlation amplitude in the x-y 
plane normal to f Q and orthogonal to the plane of Fig. 5(b) . 

Center region Ay - ± 0.5 mm was inaccessible for measurements. 

Noncircular contours indicate asymmetric azimuthal wave propa- 
gation. Preferred direction is normal to the major axes of 
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'A. -J - 


the elliptical contours. In those directions correlation 
function half-widths are comparable to those along v^. 

t_ - 24 msec, z = 7 cm, 

a 

Fig. 8. Cross correlation signal vs. distance Ax across for different 
delay times t. In contrast to Fig. 6(a) the correlation function 
splits and moves to both sides indicating wave propagation in 
opposing directions across $ . Due to oblique propagation the 
phase velocity (Vp = 2.2 x ID 5 cm/sec) exceeds sound speed; 
different correlation lengths reflect azimuthal asymmetries. 

Fig. 9. Schematic view of an array probe for observing turbulence in 
real time and space. The block diagram Indicates the method 
for recording, witn single shot exposures, six simultaneous 
traces of density fluctuations. 


Fig. 10. Single exposure of simultaneous traces of density fluctuations 

vs. time of six probes aligned in an array across B q as sketched 
on the right. Two time markers applied to all channels provide 
for an absolute time reference and scale calibration. Pertur- 
bation labelled "a" at bottom represents a growing wave propa- 

C 

gating at sound speed (c s * 1,4 x 10 cm/sec) in the direction 
from probe #6 to #1 perpendicular to v^. Perturbation "b" is a 
decaying wave. Growth and decay scale lengths are Ar c < 5 tnm 
corresponding to approximately two wavelengths; such short 
correlation lengths are characteristic for strong turbulence. 
Perturbations : w" and ,l d" propagate with velocity v > c oblique 
to the array axis. Amplitude decay along the wave crest is smaller 
than in propagation direction, Kr 2.3 x 10“^ Torr, z = 7 cm. 


< 5 Wm* 


25SS. 


* 


1 
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Fig. ll(a)Test wave interferometer signal vs. time at a constant frequency, 
applied amplitude and probe separation ( f = 3 MHz, Az ~ 12 mm). 
Propagation is along v^. Timing of current pulse is indicated on 
top. Lower trace is the boxcar sampling pulse. Oscillations 
reflect wavenumber changes due to temperature and density changes. 
Note build-up of large amplitude waves during the current flow and 
20 dB smaller signals after current turn-off. 

(b) Interferometer trace vs! position at a fixed sampling time. 

Nearly undamped submillimeter ion acoustic waves are excited from 
a line source (0.13 mm diam, 2,5 mm length) in the direction of 
the electron drift. 

Fig. 12. Frequency spectrum of the applied signal (a) and of the received 
ion acoustic wave during beam injection (b) and after end of beam 
pulse (c) {arbitrary vertical scales; t^ is the beam pulse length). 
Wave propagation in a time-varying medium gives rise to frequency 
shifts unrelated to Doppler effects or wave-wave interactions. 

Fig. 13. Test wave dispersion id vs kj j for propagation in the direction of 

the electron drift, f . ~ 3 MHz, =135 cm"”*, kr 2.3 x ICf^ Torr, 
z = 7 cm. 

Fig. 14(a)Test wave interferometer trace vs. axial position along at a 
radial distance Ax - 2 mm from the isotropic exciter at Ax = Ay = 

Az = 0. The preferred propagation in the direction of the electron 
drift is evident. 

(b)Contours of constant phases of the test wave in the plane ax (J[t ), 
az (|1 B 0 ). Circles are maxima, crosses are minima of interferometer 
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traces. Note asymmetry in parallel wavelength with respect to 
the electron drift v^ indicating the presence of an ion drift 
parallel to v^ created by ion trapping in large amplitude waves 
along . 

(c)Contours of constant relative wave amplitude. The excentricity of 
the contours with respect to the exciter at the origin indicates 
initial amplification of the diverging wave in the half space 
centered along v^. kr 2.3 x 10" 4 Tom t Q = 24 msec, z * 7 cm. 

Fig. 15. Test wave properties in a plane orthogonal to that shown in Fig. 14 
and normal to v d )[ B . (a) Interferometer traces vs ax (J_B q ) 

at different orthogonal positions Ay (J_^ 0 )* (b) Phase fronts of 
the test Wave launched at Ax = Ay = 0. (c) Contours of constant 

relative wave amplitude. Symmetric phase fronts indicate absence 
of ion drifts; asymmetric amplitude contours, previously noted in 
Figs, 7 and 8, are possibly due to radial gradients in v^. Para- 
meters as in Fig. 14. 

Fig. 16(a)Schematic probe arrangement for exciting two oblique large ampli- 
tude test waves in order to investigate their nonlinear interaction. 

(b) Interferometer traces of the two excited modes and the 

beat mode U 3 = tu-j - ug) vs. axial position along v^ at equal 
radial spacing between the two exciters. Since the off-resonant 
beat mode is destabilized by the electron drift its amplitude is 
comparable to that of the normal modes (amplitude ratio $-| : ^ : 

“ 1:2:1). t a = 24 msec, z « 7 cm, exciter probe separation 

O Q 


Ar - 3 mm. 
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Fig. 17. Phase fronts of the two test waves and their beat mode in a plane 
defined by radial and axial coordinates, Ar, Az respectively. 

The exciters are located at az = 0, Ar-| = -2 mm, Ar 2 = +4 mm. 

Points A and B are arbitrary locations at which wave vector 

■( 

matching will be demonstrated (see Fig. IB), z = 7 cm, t & = 24 msec, 
kr 2 x 10“ 4 Torr. 

Fig. 18, Examples of wave vector matching at two locations in Fig. 17. The 
three vectors ^ , It, and are determined from local wavelength 
and wave normal directions. The expected difference wave vector 
At = T<i - t 2 agrees with k 3 to within measurement accuracy. The 
phase velocity of the beat mode is, in general, smaller than sound 
speed c s and can fall into the ion distribution function (B). 
f'j Parameters as in Fig. 17. 

■•V .■ 

Fig. 19{a)Contour map of the product of the two normal mode amplitudes 
(arbitrary units}. 

(b)Contours of constant relative amplitude of the beat mode. Al- 
though the driving modes • $ 2 are largest near point B the 
beat mode amplitude maximizes elsewhere since ion Landau damping 
is strong near B. Parameters as in Fig. 17. 
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I . INTRODUCTION 

The operation of hard wire antenna systems immersed in a plasma which 
also contains an impressed magnetic field has been the subject of study 
for many years. The primary motivation has been to understand better the 
nature of signal transmission from satellites in or above the ionosphere. 

3 

For frequencies well above the electron plasma frequency Ifp e = 9 x 10 

(n e )' 5 Hz, where n g is the electron concentration in cm - ], electromagnetic 

waves propagate quite freely, since the plasma dielectric "constant," 

2 2 

e ^ 1 - (u n _/<ii , is effectively unity, and thus for ^ » w the propagation 
pe pe 

medium is approximately a vacuum. Thus, for VHF, UHF and SHF transmission, 
the predicted behavior of the antenna system may be well approximated by 
calculating the free-space radiation pattern, even though it is immersed 
in plasma. 

However, there are two conditions under which the behavior of the 
antenna system will begin to deviate from that predicted by a free-space 
theory. First, even at frequencies well above the plasma frequency, w , 
the nature of antenna behavior can be changed from free-space behavior by 
applying very hard drive (i.e., high input power) to the antenna. This 
can lead to such strong near fields that nonlinear plasma responses occur 
in the vicinity of the antenna. Examples of such apparent effects have 
been seen in data from sounder antenna experiments aboard the Alouette 
and the ISIS spacecraft (cf. bibliography section). This effect at high 
frequencies is of no concern to us in the present report. 

The second condition under which antenna behavior becomes different 
from free-space behavior is one of operation at frequencies below the 
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electron plasma frequency. For the ionosphe re/magnetos phere environment, 
this means operation at VLF, ELF and ULF. Quite aside from the problems 
associated with propagation of signals at these low frequencies in the 
inhomogeneous plasma and geomagnetic field above the it .ospheric F layer, 
and the complex wave phenomena occurring when such waves impinge upon the 
earth /ionosphere cavity, and ionospheric "waveguides 11 formed in E and F 
layers (and also above the F 2 peak), there are very important and so far 
unanswered questions concerning the way in which hard wire antenna systems 
couple energy into the surrounding ionospheric plasmas at frequencies 

til < 01 . 

pe 

II. ELECTRIC DIPOLE AMTENNAS 

Attempts to calculate the radiation from "short" electric dipole 
antennas immersed in a magnetoplasma and driven at frequencies u < 

(VLF and below) are quite numerous. They are listed under "Electric 
Dipoles" in the bibliography and reference section of this Appendix. 

The fundamental problem of calculating the radiation from an electric 
dipole immersed in a magnetoplasma is one of solving Maxwell '.s equations 
coupled to -the self-consistent field-particle equations describing the 
plasma. For very small applied antenna voltages, the plasma dynamics 
may be linearized, and the plasma reaction currents computed self- 
consistently in terms of the voltage on the antenna. Thus, the "short" 
dipole can be idealized as a point dipole of moment p Q = V Q £, where 
is the antenna element-to-plasma capacitance for a sheath of dimensions 
^ several plasma Debye lengths and £ is the dipole element's length. 

In actual practice, it is at this point that the first difficulty 
with the antenna theory arises. Even for small signals, the sheath around 
such an antenna will exhibit a complex impedance. For larger signals. 
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that Is for antenna voltages larger (or very much larger) than the mean 
thermal energy (kT & ) of an electron divided by its charge (i.e., for 
V Q > or » kT e /e) , nonlinear effects appear. The usual sheath theory, and 
the concept of linear complex impedance are no longer valid. The coupling 
of energy fed to the antenna terminals to energy actually propagated away 
from the antenna in VLF, ELF or ULF modes is quite unknown. To the present 
date, no acceptable, or even belivable theory is available to describe 
this highly nonlinear problem. 

There are rm (repeat, no_) experiments available, either in the 
laboratory or on spacecraft, dealing with high-power, VLF to ULF radiation 
into a magnetoplsama. The only experimental work available so far is that 
involying sounder antennas aboard the Alouette and ISIS satellites [see 
bibliography and reference section^, and an attempt to measure input 
impedance characteristics of a linear electric dipole on the USAF satellite 
OV 1-21 by the Aerospace Corporation group. 

Other attempts to measure impedances of various configurations of 
small dipole-like receiving antennas were made during a NASA Javelin 8.45 
rocket flight by Gurnett and co-workers at the University of Iowa. 

In no case did any of the "resonances" or other peculiar phenomena 
predicted by the available theories of dipole antennas ever occur in the 
data, even in the small signal regime. In no case did the experiment 
attempted involve any measurement of the radiation from the antenna system. 

In addition, it appears to us that no meaningful laboratory experiment 
has ever been performed to measure radiation from a dipole in a magneto- 
plasma at VLF and below. There are two laboratory experiments having some 
bearing on the problem, but they both have severe drawbacks that do not, 
in our opinion, allow extrapolation to the case of antennas in the 
ionosphere/magnetosphere. The first experiment is that of Fisher and 
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Gould ( Reference. B-l ) . 

Contrary to a commonly held misconception, the Fisher-Gould theory 
and experiment neither calculated nor measured the far-field radiation 
from a short dipole. Their theory is based on a near-zone field approxi- 
mation and involves a quasi-static (E = - v<i>) assumption. Far-zone fields 
are known to have a primary dependence on the V x E properties of Maxwell's 
equations. This same comment applies to the calculations by Kuehl 
(Reference B-2 and B-3) on the resonance cone interference phenomena found 
in oscillating electric doublets immersed in cold magnetoplasma. 

The most complete attempt to calculate far-field radiation from an elec- 
tric dipole antenna, keeping terms involving both v x E and 7 * E, has been 
carried out by Wang and Bell (References B-4 and B-5). The major drawbacks 
to these calculations are the nature of assumed plasma response (linearized 
cold^magneto-fluid) and the assumed current distribution on the antenna it- 
self (triangular). Neither assumption is generally correct, even under 
small applied antenna potentials, although the linearized fluid approximation 
is, perhaps, justifiable at small signal levels. However, the assumption 
of an effective radiating current distribution in the form of a triangular 
ramp on the two dipole elements may not so easily be justified. Statements 
to the effect that far fields are insensitive to the exact shape factor of 
the "small 11 dipole current distribution are simply extrapolations of results 
obtained in free-space calculations; there is no proof that such statements 
are correct for excitation by sources immersed in a magnetoplasma, especially 
for strong sources. 

In fact, the conditions for strongly coupling antenna energy into some 
desired mode (such as the VLF Whistler) of the magnetoplasma will be a sen- 
sitive function of the geometry of the source current distribution. Since 
no spacecraft-borne experiment has yet been performed to look at radiation 
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from either electric dipoles or magnetic loops, one can only speculate con- 
cerning the actual efficiency for radiation at VLF by such systems. It is 
our opinion that laboratory experiments are now feasible and should be per- 
formed. (See Reference B-7) 

III. MAGNETIC LOOP ANTENNAS 

For several reasons which will become clear subsequently, the magnetic 
loop antenna may perhaps avoid some of the problems besetting electric dipoles. 

The basis for this statement is simply that such loop antennas are current 
sources, rather than voltage sources, at VLF and below. The obnoxious effects 
of electrostatic (or quasi-static) sheaths and particle currents produced by 
sheath disturbances due to changing drive voltages applied to the antenna 
may be alleviated by the rather simple technical device of electrostatically 
shielding the loop conducting the current. At low enough frequencies the 
■ shield is relatively transparent to slowly varying magnetic fields, while 

preventing sheath conduction currents from becoming directly coupled to the 
antenna proper. Of course, currents may be induced in the electrostatic 
shield, and these will necessarily enter into the net driving current from 
which radiation is calculated. 

There exist, in the literature, two publications in which Wang and Bell 
(References B-6 and B-7) have treated the problem of radiation of VLF and 
ELF waves by a loop antenna immersed in a cold magnetoplasma. The basic 
assumptions employed by these authors are similar to those invoked in the 
analysis of the electric dipole radiator discussed in the previous section. 

In effect, the response of the magnetoplasma was assumed to be described 
by linearized raagneto-i.onic equations, while the loop current strength was 
given. Thus, any plasma reaction back on the loop current other than that 
} due to linear dynamical response was ignored to first approximation. 


The Wang and Bell solutions to this problem are in the form of input 
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impedence of an aribtrary orientation with respect to magnetic field 
direction of the loop's axis of symmetry. They found that the loop ori- 
entation with respect to magnetic field did not have a very strong effect 
on the input impedance at VLF and ELF. 

IV. SOME THEORETICAL ASPECTS QF ANTENNA-PLASMA COUPLING 

Most of the past theory of operation of an antenna immersed in a mag- 
netized plasma and driven at frequencies well below the electron plasma 
frequency has been based on notions that simple Debye Sheaths form, and 
plasma currents are linearized responses to the antenna's applied field. 

Although under the present contract, both time and resources prevented 
a detailed theoretical treatment of the antenna- plasma coupling phenomenon, 
it is still pertinent to outline the elements of such a theory and to point 
out the novel elements of our approach to such a theory, an approach sig- 
nificantly different from previous ones. 

A. Electric Dipole Antenna 

Consider an electric dipole antenna immersed in a magnetized plasma 
of electron and ion temperatures T T . , charge density n, and magnetic 
field B Q can be assumed. For purposes of simplicity, consider the dipole 
gap to have dimensions small compared to the length of the wires. 

One notes that a planar sheet of field lines of B intersect the 

_o 

conducting antenna elements. In the absence of any voltages applied to 
the antenna, its conducting elements will float to some small potential, 

V Q < kT e /e, ( i . s . , determined by the plasma potential), so that equal ther- 
mal currents of ions and electrons are collected by the antenna surface. 

Note that these currents are collected primarily along lines of B Q repre- 
senting guiding centers within one Larmor radius of the antenna. 


V 
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Now, consider a "d.c." operation of the dipole, in which one element 
is driven at -V-j and the other at +V^ relative to their common ground. If 
V-j exceeds the order of , then the antenna element at -V-j repels electrons, 
while that at *V- ( attracts electrons, with reverse effects on the plasma 
ions . 

Under this ''d.c. 11 mode, no change in plasma potential can occur in the 
plasma on field lines intersecting the half of the antenna at , because 
the ion current remains saturated. However, the plasma on field lines 
intersecting the half of the antenna a +V-j , electrons can be drawn from the 
plasma to the antenna, which leaves a negative charge deficiency along those 
field lines, and the plasma potential there assumes a larger positive value 
than it had previously. 

For this reason, a strong electric field will develop along an interface 
between the regions defined by field lines intersecting the positively biased 
antenna element, and regions of non-intersecting field lines. This electric 
field is perpendicular to local field lines. 

The electrons which participate in changing this plasma potential are 
predominantly fast, being the electron population with speeds on the order 
of or exceeding the mean thermal speed. Thus, a rather siziable volume 
of plasma can have its potential increased in this manner. 

If the potential on the dipole is sin mt, instead of d.c., then the 
plasma potential can be changed in volumes by each half of the antenna on 
alternate cycles, and such disturbances will lead to possible excitation 
of plasma modes, depending upon the applied frequency. 

This effect has been observed in laboratory plasmas at TRW, and re- 
quires further investigation. It appears to represent a mechanism whereby 
excitation of a plasma at frequencies well below the electron plasma 
frequency can be effected. 
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A general model based on this concept can be constructed and examined 
theoretically and investigated experimentally. Because the electric fields 
produced are predominatly normal to local B Q > it appears certain that very 
effective excitation of Bernstein and cyclotron waves will occur under 
appropriate conditions on the frequency. 


I 
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WHISTLER SHOCK RELAXATION 


Background 

One of the modes of electromagnetic wave propagation in the solar wind 
is the whistler mode familiar in the magnetosphere and the laboratory. Pre- 
sumably, any source of wave generation in the solar wind should produce 
whistlers, which will travel at phase and group velocities dependent on fre- 
quency, on the local plasma parameters, and on the direction of propagation with 
respect to the interplanetary magnetic field {the IMF). Any shock wave, In- 
cluding the earth's bow shock, is such a source, and the detection and behavior 
of whistlers in the solar wind near the bow shocks of Earth and Mercury has 
been the subject of reports by Fairfield (197^) and Fairfield and Behannon ' 

0 976). 

In general, waves cannot propagate upstream from the bow shock Into the 
oncoming solar wind. Indeed, the shock forms because the wind flows "super- 
magnetosonically." However, whistler, or fast, waves constitute the one mode 
that includes zero frequency in its range and that, for some frequencies and di- 
rections, can propagate upstream under fairly common parameter conditions. 
Dispersed whistler phase and group velocities ar highest, regardless of 
direction, at frequencies somewhat below the elect qyrof requency f . 

These are the frequencies at which wave speeds are most likely to exceed the 
solar wind velocity V sw , with waves running away from the shock, as reported 
by the authors cited above. At lower frequency, the phase velocity may match 
Vgy, with the result that the perpendicular-shock jump profile is strongly 
modified by the addition of whistler standing waves ahead of It In the solar 
wind, at least in laminar structures, when the IMF is about 10°-40° from the 
shock "front. 11 Such standing wave profiles have been observed in the laboratory 
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( Robson , 1969) and in space ( Greenstadt et al., 1975; Fairfield and Feldman , 

1 975) » The whistler mode, then, is an important component, both theoretically 
and experimentally, in shock structure. 

Certainly the most profound significance of the whistler mode for shock 
structure must occur when waves with the widest possible band of frequencies, 
especially including those low enough to approach the Ion cyclotron frequency, 
are capable of overmatching the solar wind velocity and carrying a significant 
portion of the shock’s energy upstream. This situation probably contributes 
to the total breakdown of shock discipline when the IMF, B, Is less than about 
50° from the local shock normal ( Greenstadt et al . , 1976). The properties of 
whistlers that lead to such shock "relaxation" is the subject of the study 
reported in this section. 

Relaxation of a quas i -paral 1 el plasma shock is a process of considerable 
complexity involving counterstreaming particles, local electric and electro- 
magnetic wave excitation, wave amplification, reduced stream thermal ization, 
upstream plasma modification, and energy escape by wave propagation. The 
part played by simple wave escape through the whistler mode has never been 
Investigated thoroughly, particularly for typical solar wind plasmas, arbi- 
trary angles of the IMF, and realistic curved shock profiles. 

Objective 

The objective of this study is to characterize fully the behavior of 
the whistler mode for arbitrary angles of propagation to B, under representa- 
tive solar wind conditions. It is hoped that in this way the contribution of 
whistlers to quasi-parallel structure can eventually be assessed and des- 


cribed. 
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Approach 

The Initial attack on the problem outlined above proceeds from a simpli- 
fied basis with the intention of restoring neglected complications later. 

The sections of the whistler dispersion curves where dissipation of shock 
energy must take place first are those where phase and group velocities are 
highest, i.e., where f . « f < f . In this range of frequencies the 
Appleton-Hartree (A-H) approximation applies for cold plasma, and th/s defines 
the beginning approach. 

The A-H dispersion relation for the fast mode ( St?x , 1962, p kO) gives 
the phase velocity expression 



where u, = u,k, 9 is the angle between propagation vector k and B, and f 

**(p * - ~ ce 

and f are the electron cyclotron and electron plasma frequencies. The group 

velocity u of the waves differs in general in both speed and direction from 
9 

□ , . if we take 
~(j) 

u = u cos 0B+u sin0Bi 

~g g g- g s\L 


where 9 is the angle between u and B, and Bi is a vector perpendicular to 
g ~9 ~ ~± 

d to 

B in the plane of u and B, the selected approximation gives, from u = 


u z = 


J d> 


cos^S 


3 cos z 0 


- 8 j— cos 0 + k (p— \ + 1 1 
ce V ce / J 


tan 0 

g 


(cos 0-2 



cos 0-2 



1 

cos 0 


and 
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Obviously, even this simplified cold plasma approximation connects u 
0^ through expressions whose properties are less than self-evident. 


9 


and 


The plan is to exploit the above approximation to study the properties 
of whistler propagation in realistic solar wind models surrounding a model 
bow shock to see if any pattern of wave escape can be linked to experimental 
results describing the development of quas i -para 1 lei structure. 


Status 


In addition to derivation of the auove mathematical model, a program 
has been written to compute and plot various versions of the elaborated re- 
lationships, and preliminary graphs have been constructed for a few parameter 
conditions, principally to test the program. Figure 1 Is a sample polar plot 

of phase ano group velocities vs^ angle for one frequency and one set of 

"3 

typical (colo) solar wind quantities, B = 5y, N = 7 cm , which determine 
f ce /f . The figure tells us that, for the selected quantities, the phase 
velocity would exceed a nominal solar wind velocity of 400 km/sec (dashed 
curve) at all angles 0 < 52°, while the corresponding group velocity would 
be more than double the solar wind speed, but only for angles 6^ £ 15°. 

In application, the actual relationship of B to the local shock normal 
and the appropriate component of solar wind velocity must be taken into 
account to determine the pattern of whistler wave escape from the shock. At 
the time o' 7 writing, the study has not been carried beyond production of a 
few test plots such as that In Figure 1. Eventually, it is hoped that a 
comprehensive picture of wave front and bow shock geometry will be developed. 
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. Polar plot of whistler phase and group velocities 
for one set of solar wind parameters. 
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